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1. Objectives 

The purpose of this meeting is to discuss the current status and future development of 
turbulence modeling in computational fluid dynamics for aerospace propulsion systems. 
Various turbulence models have been developed and applied to different turbulent flows over 
the past several decades and it is becoming more and more urgent to assess their perfor- 
mance in various complex situations. In order to help users in selecting and implementing 
appropriate models in their engineering calculations, it is important to identify the capabil- 
ities as well as the deficiencies of th ese models. This also benefits turbulence modelers by 
permitting them to further improve upon the existing models. 

This workshop is designed for exchanging ideas and enhancing collaboration between differ- 
ent groups in the Lewis community who are using turbulence models in propulsion related 
CFD. In this respect this workshop will help the Lewis goal of excelling in propulsion related 
research. 

2. Format 

• This meeting has seven sessions for presentations and one panel discussion over a period of 
2 days. 

• Each presentation session is assigned to one or two branches (or groups) to present their 
turbulence related research work. Each group should address at least the following points: 
current status of turbulence model applications and developments in the research; progress 
and existing problems; requests about turbulence modeling. 

• Each speaker will be given 18 minutes for the presentation which will be followed by 2 
minutes of questions/answers. 

• The panel discussion session is designed for organizing committee members to answer man- 
agement and technical questions from the audience and to make concluding remarks. 

• All the talks will be collected and printed in a proceedings of this workshop. 

3. Organizing Committee 


L.A. Povinelli (Chairman) 
J.M. Abbott 
P.G. Batterton 
R.E. Gaugler 


T. Keith 
R. Mankbadi 
E.J. Mularz 
D.R. Reddy 


R.J. Shaw 
R.J. Simoneau 
R.M. Stubbs 
A. Shabbir 
T.-H.Shih 


1 



Workshop on Computational Turbulence Modeling 
OAI Auditorium 
September 15. 1993 


8:40 

/- 9:00 
2 - 9=30 
O 9:50 
J 10:10 
p 10:30 

O 10:50 
o 11:10 

q 11:20 
12:00 



&3:20 
tf 3:30 

C O 3:50 
A; 4:00 

s/ 

P 4:30 


Welcome by L. Povinelli 

Turbulence Modeling Activities at CMOTT 
Chairman: T.-H. Shih 
CMOTT Research Activities by T.-H. Shih 
Research on two equation models by Z. Yang 
A Reynolds stress algebraic equation model by J. Zhu 
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PDF models by A. Hsu (CFD branch) and A. Norris 
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Two-equation models for turbine blade heat transfer by A. Ameri 
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Introduction: Development of New Flow Reactive Code (ALLSPD) by E. Mularz 

A coupled implicit solution method for turbulent spray combustion in propulsion systems 
by K.-H. Chen and J.-S. Shuen 

Research Activities by R.M. Stubbs 

On turbulent transport of chemical species in compressible reacting flows and Unsteady 
transitional flows over forced oscillatory surfaces by S.-W. Kim 

On the accuracy of compact differencing schemes for DNS by S.-T. Yu 
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Research Activities at the Center for / £ &¥& "3 
Modeling of Turbulence and Transition /°* 

Tsan-Hsing Shih 

1. Introduction 

2. General developments 

2.1 Turbulent constitutive relations 

2.2 Mechanical and scalar dissipation equations 

2.3 Eddy viscosity transport equation 

3. One-point closure schemes 

3.1 One-equation eddy viscosity transport equation model 

3.2 Galilean and tensorial invariant realizable k-e model 

3.3 Reynolds stress algebraic equation model 

3.4 Scalar flux algebraic equation model 

3.5 Reynolds stress transport equation model 

3.6 Non-equilibrium multiple- scale model 

3.7 Bypass transition model 

3.8 Joint scalar PDF model 

4. RNG and DIA 

5. Numerical simulation 

Abstract 

The main research activities at the Center for Modeling of Turbulence and Tran- 
sition (CMOTT) are described. The research objective of CMOTT is to improve 
and/or develop turbulence and transition models for propulsion systems. The flows 
of interest in propulsion systems can be both compressible and incompressible, 
three dimensional, bounded by complex wall geometries, chemically reacting, and 
involve “bypass” transition. The most relevant turbulence and transition models for 
the above flows are one- and two-equation eddy viscosity models, Reynolds stress 
algebraic- and transport-equation models, pdf models, and multiple-scale models. 

All these models are classified as one-point closure schemes since only one-point (in 
time and space) turbulent correlations, such as second mo ment s (Reynolds stresses 
and turbulent heat fluxes) and third moments (u t UjU k , u;0 2 ), are involved. In 
computational fluid dynamics, all turbulent quantities are one-point correlations. 
Therefore, the study of one-point turbulent closure schemes is the focus of our 
turbulence research. However, other research, such as the renormalization group 
theory, the direct interaction approximation method and numerical simulations are 
also pursued to support the development of turbulence modeling. 
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1. Introduction 

The center for modeling of turbulence and transition was established as a special 
focus group within the Institute for Computational Mechanics in Propulsion at 
NASA Lewis Research Center in 1990. Its objective is to improve and/ or develop 
turbulence and transition models for computational fluid dynamics (CFD) applied 
in propulsion systems. With the advance of computer technology and algorithms, 
accurate turbulence and transition modeling becomes the pacing item for improving 
flow calculations used in propulsion system design in all its key elements. The flows 
of interest in propulsion systems are, in general, very complex since there are wall- 
bounded three-dimensional complex geometries, chemical reactions, compressibility 
and transition, etc. In order to accurately predict these flows one must correctly 
model the turbulent stresses and scalar fluxes which are one-point (in time and 
space) turbulent correlations. For flows with finite rate chemical reactions, accurate 
modeling of the production rate of species is crucial for turbulent flow calculations. 
Based on the above considerations, turbulence modeling activities at CMOTT are 
focused on one-point closure schemes, that is, using the moment closure schemes for 
the turbulent velocity field and the joint scalar pdf method for the reacting scalar 
field. 

There are various moment closure schemes which have been developed for var- 
ious engineering applications. However, in practice, one often finds that the ex- 
isting models need to be improved and/or re-developed in order to reasonably 
simulate complex flow structures appearing in propulsion systems. For this pur- 
pose, CMOTT devotes itself to improving and/or re-developing these moment clo- 
sure schemes which include eddy viscosity (one- and two-equation) models, second 
moment algebraic- and transport-equation models, non-equilibrium multiple-scale 
models, and bypass transition models. In addition, other studies supporting the 
development of one-point closure schemes have been also carried out (for example, 
studies on renormalization group theory (RNG), direct interaction approximation 
(DIA), direct numerical simulation (DNS) and large eddy simulation (LES)). 

In this report, we first describe the general development of turbulent constitutive 
relations, turbulent mechanical and thermal dissipation and a new eddy viscosity 
equation. Second, we describe the detailed developments on each moment closure 
scheme and the pdf method. Then the RNG and DIA methods and finally, the 
numerical simulation of particular turbulence phenomena, such as rotation and 
bypass transition, etc., are considered. 

Each research subject is the joint project of several CMOTT researchers and 
visitors. In describing research activities, the names of involved researchers will be 
mentioned for reference. 
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2. General Developments 

2.1 Turbulent Constitutive Relations 


Reynolds stress 

Using the invariant theory in continuum mechanics and Generalized Cayley- 
Hamilton formulas for tensor products, a turbulent constitutive relation (or a gen- 
eral turbulence model) for any turbulent correlations can be obtained, in principle. 
Therefore, this theory provides an avenue to develop better turbulence models than 
those existing. For example, a commonly used constitutive relation for Reynolds 
stresses uiu] (in terms of the mean deformation rate tensor TJ^j and the turbulent 
velocity and length scales characterized by the turbulent kinetic energy k and its 
dissipation rate e) is 


2 

—U{Uj = Cp {U{j + Uj y i) — -k6ij 
£ o 


The effective eddy viscosity vt defined as 


-UiUj k 2 

VT = u-Tu-i= C ^ for 


( 2 . 1 . 1 ) 


(2.1.2) 


is isotropic since vt is a scalar quantity. However, the invariant theory enables us 
to formulate the following general model (Shih and Lumley 1 , Johansson 2 ): 

2 K 2 2 

UiUj — kS 2 j T 2a 2 ~ "I" 

+ 2a i ^-{Vfj + Ul-^a ,«„) 

+ 2 a 6-~2~(Ui,kUj,k — rll2 Sij) 

£ o 

K z 1 

2^7 ~^2~ {Uk,iUk,j — ~n2 Sij) 

+ 2a s tA {Ui ^ jk + V f k Ujk _ |n 3 ^.) 


+ 2a l0 ^.(u k<i ul d + UtjUl - - 3 il 3 %) 

K 5 1 

+ 2a 12 -(Ul k Ul k --U 4 6 tJ ) 


K 5 1 

AUlJJ 2 
£ q 


(2.1.3) 


K 5 2 

+ 2au—(Ui'kUi t kUij + U^U^Vl - 

K 6 9 

+ 2016— (Uij.Vb.itij + u^u^ul - |n 6 « 0 .) 

+ 2o 18 ^(ci i , i n,, t pf iin ^ ra + Vu.Uij.vljJl, - ?n 7 ^) 


4 


Tsan-Hsing Shih 


where 


III — , 

n s = Ui+u^ul 


II2 — Ui^kUi'k i 


n 3 


n 6 = Ui,kUl k ul 


Ui, k ult , n 4 = uyjh ■ 
n 7 = u^v^ulJJln 

(2.1.4) 


From Eq.(2.1.3), the effective eddy viscosity 


( u T)ij - 


— UiUj 

Uij + Uj'i 


(2.1.5) 


is no longer a scalar and, hence, is an anisotropic eddy viscosity. It is noticed that 
the first two terms on the right hand side of Eq.(2.1.3) represent the standard k-e 
eddy viscosity model (2.1.1) and that the first five terms of Eq.(2.1.3) are of the same: 
form as the models derived from both the two-scale DIA approach (Yoshizawa 3 ) and 
the RNG method (Rubinstein and Barton 4 ). 

Eq.(2.1.3) is a general model for uiu ] . It contains 11 undetermined coefficients 
which are, in general, scalar functions of various invariants of the tensors in ques- 
tion, such as SijSij (strain rate) and Q,ijQ.ij (rotation rate) which are (H 2 "I - Hi)/^ 
and (II 2 - IIi)/2 respectively. The detailed forms of these scalar functions must 
be determined by other model constraints, for example, realizability, and by exper- 
imental data. Eq.(2.1.3) contains 12 terms; however, its quadratic tensorial form 
may be sufficient for practical applications. We will see later in section 3.3 that the 
constitutive relation (2.1.3) has a significant impact on the development of Reynolds 
stress algebraic equation models. 


Turbulent scalar flux 8u{ 

We assume the following functional form: 


Qui = Fi(Uij,T t i,k,£,9 2 ,£q) 


(2.1.6) 


where 0 2 is the variance of a fluctuating scalar and £0 is its dissipation rate. 
Eq.(2.1.6) indicates that the scalar flux depends on not only the mean scalar gradi- 
ent T j, but also the mean velocity gradient Uij and the scales of both velocity and 
scalar fluctuations characterized by k, e, 0 2 , £ 0 - 

Applying the invariant theory, we may obtain the following general constitutive 
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relation for Qui'. 


Wi = aifc(--) 1 / 2 r,i + 

v e ee ' e e e$ 

+ — r (- — ) 1/,2 (tt4^i,fc^fcj + 0'$Uj'kUk,i + CL6Ui t kUj,k + a7Uk,iUk,j )Tj 
e l see 

+ ^ (- — ) 1/2 (asl7i ) jfcl7| fc + ag^,* +o l0 l7 fc>i Oj ii + an^Dikj )Tj 


e* £ £<? 


+ auUitUwUtj + ai$U j,kUi t kUi,i )T,j 



* Se ° (2.1.7) 

The coefficients ai - ais are, in general, functions of the time scale ratio -/ — 
and the other invariants formed by the tensors in question, for example, T k T k , 
TiUijTj, etc.. Again, Eq.(2.1.7) implies that the effective eddy diffusivity 


(7t)» = 



is not isotropic. It is noticed that the first term on the right hand side of Eq.(2.1.7) 
is the standard eddy diffusion model, and the models derived from the two-scale 
DIA (Yoshizawa 5 ) and the RNG method (Rubinstein and Barton®) are similar to 
the first two terms of Eq.(2.1.7). In practice, a form containing the first two terms 
on the right hand side of Eq.(2.1.7) may suffice. Further development of this model 
for turbulent heat transfer is described in Section 3.4. 

The Researchers involved with the subject in this section are T.-H. Shih, J. Zhu, 
A. Shabbir, J.L. Lumleytand A. Johansson.^ 

2.2 Mechanical and Scalar Dissipation Equation 


Mechanical dissipation e 

In turbulence modeling, we often need turbulent characteristic velocity and length 
scales. While the turbulent kinetic energy k is used to characterize the velocity scale, 
the mechanical dissipation rate e and the scalar dissipation rate ee are used to char- 
acterize the length scales for mechanical and scalar fields, respectively. Comparing 
with the turbulent kinetic energy equation, the exact dissipation rate equation is 


f Professor, Cornell University, Ithaca, NY 
| Professor, Royal Institute of Technology, Stockholm, Sweden 
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very complicated. In this equation, all the terms which represent important tur- 
bulence physics (for example, turbulent diffusion, generation and destruction) are 
unknown and are of complex forms that are all related to small scales of turbulence. 
Therefore, in the literature, the exact dissipation equation is not considered as a 
useful equation to work with. Instead, one creates a model equation by assuming an 
analogy to the turbulent kinetic energy equation, i.e., one assumes that the model 
dissipation rate equation also has generation and destruction terms which are as- 
sumed to be proportional respectively to the production and dissipation terms in 
the turbulent kinetic energy equation over the period of large eddy turn-over time 
characterized by kje. The resulting model dissipation rate equation is written as 

£,t + Ui£j =V£,U — (€Ui)g 

£ £? 

Cel~Ui'UjUi > j 

Recently, Lumley 7 proposed a dissipation rate equation based on the concept of 
spectral energy transfer caused by interactions between eddies of different sizes. 
This model equation mimics the physics of statistical energy transfer from large 
eddies to small eddies and is of a different form than equation (2.2.1). 

In this study, we explore another rational way to obtain the model dissipation 
rate equation which contains certain important physics and hope it will work better 
than the existing one. The idea is that first, there is a relationship between the 
dissipation rate e and the mean-square vorticity fluctuation uiul at high Reynolds 
numbers or in homogeneous turbulence: 


( 2 . 2 . 1 ) 


£ = VUiUi 

and second, all the terms appearing in the UiuTi equation have more clear physical 
meanings than that in the £ equation so that the uJiul equation is easier to model. 
Once the UJui equation is modeled, a model dissipation rate equation will be readily 
obtained. 

The exact equation for uTgUl is 


2~ + u j{—2~),j = K-j -Ud ~ 2 {ujUiWi),j + UiUijSlj 

- UjuJiQij + l uiUjUij + cxJiUJjUij - vtOijuJi.j 


( 2 . 2 . 2 ) 


where Ui and U l are the fluctuating and mean velocities, and u)i and are the 
fluctuating and mean vorticity which are defined by 


k^k,J f “ i — £ijkUk,j (2.2.3) 

Tennekes and Lumley 8 clearly described the physical meaning of each term in equa- 
tion (2.2.2). Order of magnitude analysis shows that the first, third, fourth and fifth 
terms on the right hand side of Eq.(2.2.2) become small compared with all other 
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terms in the equation as the turbulent Reynolds number increases. The sixth and 
seventh terms are the production due to fluctuating vortex stretching and the dissi- 
pation due to the viscosity of the fluid. As the turbulent Reynolds number increases 
these last two terms become dominant and the balance between them determines the 
evolution of vorticity fluctuations. Neglecting terms uiuijQj, ul HJjUij 

and the evolution at large Reynolds number will be described by 

the following equation, 


/ \ yt i \ / \ , 

( 2~ X* + ^i( 2 ~~ )>i “ WiWjUij — utOijUij 


To model uiiUjUij — vuijaJij, let us first estimate UiUjUij. 
anisotropic tensor 


1 

~ ul O 3 


(2.2.4) 
we define an 

(2.2.5) 


then UiUJjUij can be written as 


UiUjUij = btyluij 


( 2 . 2 . 6 ) 


We expect that the vortex stretching tends to align vortex lines with the strain rate 
so that the anisotropy bfj would be proportional to the strain rate s t j, i. e., 

bfj oc — , where s = (2sijSij) 1 ^ 2 , Sij = (uij + Uj ^)/ 2 (2.2.7) 

s 


This leads to the following model: 


OJiUJjUi'j OC 0^ \j 


( 2 . 2 . 8 ) 


where we have assumed that and (2sijSij) 1 ^ 2 are well correlated. 

Using the relation, u)i = eijkUkj, it is not difficult to show that at large turbulent 
Reynolds number, 

uJiuTi ~ 2 SijsTj (2.2.9) 

and Eq. (2.2.8) can be also written as 



( 2 . 2 . 10 ) 


Equation (2.2.10) indicates that this term is of the order (u 3 /l 2 )R^ 2 as it should 
be. On the other hand, from eq. (2.2.4) the term uiaJjUiJ — must be of the 

order ( u z /l 3 )R t which is the order of magnitude of all the other terms in Eq.(2.2.4), 
therefore the term —vZJ-JuTj-j must cancel the term (2.2.10) or (2.2.8) such that the 
difference of these two terms is smaller than the term (2.2.10) or (2.2.8) by an order 
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of . This suggests that the combination cj z iOjU l j - vutijWij can be modeled by 
the following two terms: 


“I “i 


k + 


a AS 


( 2 . 2 . 11 ) 


because the ratio of k/ v to and the ratio of s to S are of order where 

k « u 2 is the turbulent kinetic energy and S is the mean strain rate (2 SijSij) 1 / 2 . 
Equation (2.2.11) does give the right order of magnitude for 
Therefore, the dynamical equation for fluctuating vorticity (2.2.4) at large Reynolds 
number can be modeled as 


(^),t + = ~(t Tjuimh + c ul uA k s - c* 


■ OJiUJi 


<*>l <d 


k + 
1/ ' 


( 2 . 2 . 12 ) 


Using e = vW h we readily obtain the following model dissipation rate equation, 

e,t + Uj£,j = —(u]e),j + C^iS e - C w2 ^ ( 2 . 2 . 13 ) 

where C w \ and C w2 are the model coefficients which are expected to be constant at 
large Reynolds number. 

It should be noticed that Eq.(2.2.13) is different from the standard £ equation 
(2.2.1) by both the generation and destruction terms. First, the Reynolds stresses 
do not appear in the generation term and the new form of the generation term is 
similar to that proposed by Lumley 7 which is based on the concept of spectral energy 
transfer. Second, the destruction term is well behaved so that equation (2.2.13) will 
not have a singularity anywhere in the flow field. We expect that equation (2.2.13) 
will be numerically much more robust than equation (2.2.1). 

Equation (2.2.13) can be applied to any level of turbulence modeling including 
second order closure models; however the turbulent transport term (euf),* needs 
to be modeled differently at different levels of turbulence modeling. In an eddy 
viscosity model, the term ( eui),i will be modeled as 

(m), i = -(— e,i),i (2-2-14) 

The coefficients C w i, C w2 , cr e and the eddy viscosity v T must be calibrated using 
experimental data (Shih et al. 9 ) 


Scalar dissipation £q 

A similar analysis leads to the following model scalar dissipation rate equation: 


£e,t + Uj£e,j = ~(uj£' e ),j + C^S £$ + C^Pr 1 ^ §^££8 - Cg 3 - — -—= (2.2.15) 


where $ = y/T^Fj and T is the mean scalar quantity, such as, the mean tempera- 
ture. Further development of heat transfer model is described in Section 3.4. 

The Researchers involved with the subject in this section are T.-H. Shih, W. Liou, 
A. Shabbir and Z. Yang. 
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2.3 Eddy Viscosity Transport Equation 

In eddy viscosity models, one accepts the following simple constitutive relation 


UiUj — 2z/j - Sjj "f~ ^ kS i; j 


(2.3.1) 


and assumes that the eddy viscosity is characterized by some kind of velocity and 
length scales u' and l: 

vt ocu' l (2.3.2) 

In two-equation k-e eddy viscosity models, for example, one specifies that 


, ,i , kl 

U OC , t oc — 


(2.3.3) 


and, hence, the eddy viscosity is assumed as 


vt — cj~~ (2.3.4) 

The eddy viscosity assumption (2.3.4) is commonly adopted in two-equation models. 
Eqs.(2.3.1) and (2.3.4) together with appropriate k and e equations have been widely 
used in engineering calculations. However, for cases where the mean flow changes 
quickly or has a strong mean stream-line curvature or rotation, etc., this kind of 
model does not work very well, since the assumption (2.3.4) is too simple to account 
for the effect of the above mean flow structure on eddy viscosity. 

The main purpose of this study is to drop the assumption (2.3.4) and to derive 
an exact equation for v T based on Eq.(2.3.1) and other exact turbulence equations 
(i.e. first principles). In this way, we hope that some important turbulent physics 
can be brought into the eddy viscosity and that a physically sound turbulence eddy 
viscosity can be calculated. 

Using Eq.(2.3.1), we may write for incompressible flows 

uWj UiUj = 2 ^ k 2 , where S 2 = 25^-5^ (2.3.5) 

O 

Differentiating both sides, we obtain 



Sij D V T D „n 

S 2 Dt UiUj 2 S 2 Dt S 


(2.3.6) 


The equation for UiUj can be written as 
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where 

Dij = [vuiUj - UiUjUk],k 

Pij = -WjLkUj,k - UjUkUi^k 

„ 1 

% = — pm + p,jUi 

p 

£ij — %l / Ui j k'U'j,k 

COij — — 2£t TO fcfij7j'Ufc'Uj 2 Ejmied-rnUkUi 

Inserting Eq.(2.3.7) into Eq.(2.3.6), we obtain an exact transport equation for eddy 
viscosity 

^t = -§f(A; + Pa + n,y - ea + Coy) - ( 2 - 3 * * - 8 ) 

In this equation, all the important turbulence physics in the Reynolds stress equa- 
tion, such as Reynolds stress diffusion term Dij , production term Pij, pressure- 
velocity gradient correlation term 11^ and dissipation tensor £ij, are involved. 
Comparing with the standard eddy viscosity assumption (2.3.4), this exact eddy 
viscosity equation (2.3.8) contains very rich turbulence physics. This equation also 
im plies that a second order closure model will naturally lead to a corresponding 
eddy viscosity model. 

Now, as an example, we use Launder Reece and Rodi’s 10 model and a gradient 
transport model for the triple velocity correlation (—TfiUjUk = ^ u i u j^ t0 derive 
a model equation for vt- The resulting equation is 


3 VT . VT,kS% 


D. Vt = [(„ + + (*+ S-) 


+ k — C’i - V'i + 2(C2 — 
15 Jv 


S 2 

SikSkj Sji 

52 


vt 2 vxSij Sij^kk 

~o r) 52 


vt 


D 


(2.3.9) 


2S 2 Dt 


Note that the Coriolis terms do not explicitly appear in this equation; however 
the rotation effect on v T could be carried over through the mean flow field. In 
addition, we also note that there are no extra model coefficients introduced in 
Eq.(2.3.9). All model coefficients ( a , Ci and C 2 ) are brought in from the second 
order closure model. The values of these model coefficients may need adjustment in 
model applications. Note that Eq.(2.3.9) is not a self-consistent equation since the 
turbulent kinetic energy k and its dissipation rate £ are also involved. Eq.(2.3.9) 
together with k-£ transport equations will provide a new three-equation model which 
may better represent the effect of mean flow structure as well as mean flow history 
on the eddy viscosity. 

The Researchers involved with the subject in this section are T.-H. Shih, Z Yang, 
and W. Liou. 


3. One-Point Closure Schemes 

In this section, we describe the developments on each of the moment closure 

scheme and the pdf method which are of concern at CMOTT. The first two sections 
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3.1 and 3.2 describe the one- and two-equation isotropic eddy viscosity models. Sec- 
tions 3.3 and 3.4 describe the new developments on Reynolds stress and scalar flux 
algebraic equation models. Section 3.5 assesses Reynolds stress transport equation 
models. Section 3.6 describes a multiple-scale model for non-equilibrium turbulence. 
Section 3.7 is about transition models. Finally, in Section 3.8 the pdf method for 
turbulent chemical reaction is described. 

3.1 One-equation eddy viscosity model 

Recently developed one-equation eddy viscosity models are either based on the 
assumption (Baldwin and Barth 11 ): 


i >t — (3.1.1) 

or created according to computational experience (Spalart and Allmaras 12 ). Both of 
them are successful in some flow calculations. This scheme is quite attractive in CFD 
because one only needs to solve one scalar vt equation without bothering about 
other turbulence quantities. However, comparing with k-e two equation models, the 
above mentioned one-equation v? models do not contain any more turbulent physics. 
In fact, Baldwin and Barth’s model is, basically, a change of dependent variable 
based on Eq.(3.1.1) plus some extra approximations. Therefore, in principle, we 
should not expect any superior performance over two-equation models. However, if 
we do not use the assumption (3.1.1), there is the possibility to improve and extend 
the capability of one-equation eddy viscosity models. 

The objective of this study at CMOTT is to derive a physically sound eddy 
viscosity equation which contains rich turbulent physics and accounts for various 
effects from mean flow structures. 

Note that in Section 2.3 we have already derived an exact equation for the eddy 
viscosity (2.3.8) and also a model equation (2.3.9) which is based on the Reynolds 
stress transport equation model of Launder, Reece and Rodi (LRR). All turbulent 
physics contained in the Reynolds stress equation can be brought into the eddy 
viscosity equation. Therefore, in principle, the transport equation (2.3.9) should be 
better than existing one-equation models based on Eq.(3.1.1). However, Eq.(2.3.9) 
is not self-consistent because k and £ are also involved. To make Eq. (2.3.9) self- 
consistent, we must model k and k/s in terms of vt and 5. In most shear flows, 
the energy-containing eddy turn-over time kje is of the same order as the mean 
flow time scale S -1 , so that e/k cc S is a reasonable model. In addition, a crude 
dimensional analysis gives k oc vtS and this is, of course, reasonable only for shear 
flows. After the above considerations, the resulting self-consistent one-equation 
model is: 


D r/ , v T\r \ i CU 
-v T = [(V + T )W,tU - — 


+ 2(C„2 — 2)1/7’ 


SikSkjS 

S 2 


Jl 


{VT),k( u T),k + C V \S Vt 

vt D S 2 
' 2 S 2 Dt 


(3.1.2) 
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where the diffusion terms from the Reynolds stress equation (2.3.7) have been ma- 
nipulated and approximated. Eq.(3.1.2) clearly exhibits the various effects of the 
mean flow on the eddy viscosity. 

The model coefficients C u \ , C u2 and o can be determined by using the experi- 
mental data of homogeneous shear flows, free shear flows and boundary layer flows 
as well as the relations in the inertial sublayer. Extensive tests of this model in 
various flows are carrying out at the CMOTT. 

The Researchers involved with the subject in this section are T.-H. Shih, W. Liou, 
Z. Yang and J. Zhu. 

3.2 Galilean and tensorial invariant realizable k-e model 

The two-equation k-e eddy viscosity model is one of the most widely used tur- 
bulence models in engineering calculations. The k-e model has versions for high 
Reynolds numbers and for low Reynolds numbers. For wall bounded turbulent 
flows, the high Reynolds number k-e model (for example, Launder and Spalding 13 ) 
must be applied together with a wall function as its boundary condition, while the 
low Reynolds number k-e model (for example, Jones and Launder 14 ) can be inte- 
grated to the wall. The high Reynolds number k-e model of Launder and Spalding 
is considered as a standard k-e model. We notice that even though the model dissi- 
pation rate equation is created by assuming an analogy with the turbulent kinetic 
energy, there was not much modification until Lumley 7 and Shih et al. 9 For near 
wall turbulence, in addition to Jones and Launder’s model, there are many other 
versions of low Reynolds number k-e models (such as Chien 15 , Shih and Lumley 16 , 
Yang and Shih 17 ) which have made better performance over Jones and Launder’s 
model. 

There are, probably, four or five issues worth mentioning about existing low 
Reynolds number k-e models: the model constants are not consistent with those 
in the high Reynolds number k-e model; the wall correction terms and damping 
functions are related to the wall distance so that models are not tensorial invariant; 
a nonrealistic dissipation rate near the wall is introduced; they are not always 
realizable since normal stress could become negative; and finally, they do not work 
very well for boundary layer flows with various pressure gradients. 

The objective of this study at CMOTT is to overcome the above mentioned 
problems. First, we propose a vorticity dynamics based dissipation rate equation 
as a part of high Reynolds number k-e base model. 9 Second, based on the invariant 
theory, inhomogeneous terms for the dissipation rate equation are proposed which 
enable the model to better respond to the change of pressure gradients (Yang and 
Shih 18 ). Third, the wall distance parameter is removed from the damping function 
so that the model is tensorially invariant (Yang and Shih 19 ). The model constants 
are consistent with those in the high Reynolds number k-e model. Finally, the 
non-negativity of normal Reynolds stresses, the realizability condition, is imposed. 

The Researchers involved with the subject in this section are Z. Yang, T.-H. Shih 
and C. Steffen. 
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3.3 Reynolds stress algebraic equation model 

All eddy viscosity models including one- and two-equation models are isotropic. 
For the flows where anisotropy is important, for example, the secondary flows driven 
by turbulent normal stresses in a square duct or curved duct, eddy viscosity mod- 
els do not produce correct flow structures. To overcome this intrinsic deficiency of 
isotropic eddy viscosity models, one proposes a Reynolds stress algebraic equation 
model which will provide an effective anisotropic eddy viscosity. The first such a 
model was proposed by Rodi 20 and it achieved some success in the prediction of 
anisotropic related flow structure. However, Rodi’s formulation is a set of algebraic 
non-linear system equations for Reynolds stresses and it often creates numerical dif- 
ficulty in obtaining a converged solution. Recently, Taulbee 21 obtained an explicit 
algebraic expression for the Reynolds stress using Pope’s 22 tensor expansion formu- 
lation and solved this numerical difficulty. However, in general, Rodi’s formulation 
assumes that the ratio U{Uj / k is constant and, of course, this is not really true for 
most turbulent flows of interest. Therefore, sometimes, this Reynolds stress alge- 
braic equation model produces even worse results than the isotropic eddy viscosity 
models for cases where eddy viscosity models are appropriate. 

Alternative ways for obtaining effective anisotropic eddy viscosity models have 
been tried by a few researchers, for example, the DIA method by Yoshizawa 3 , 
the RNG method by Rubinstein and Barton 4 and invariant theory by Shih and 
Lumley. 1 . It is interesting to point out that the RNG and DIA methods result in 
the same formulation and that this formulation is the first five terms of a general 
constitutive relation Eq.(2.1.3) except that the model coefficients are different. 

One of our goals at CMOTT is to search for an effective anisotropic eddy vis- 
cosity model for complex turbulent flows where the nonequilibrium of turbulence 
is not very severe so that the constitutive relation (2.1.3) is more or less valid. We 
have explored the potential capability of Eq.(2.1.3) and found that a truncation of 
Eq.(2.1.3) up to the quadratic terms of the mean velocity gradients is sufficient for 
various flows of interest. The model coefficients are determined such that realizabil- 
ity for the normal stresses is ensured. The detailed analysis is described by Shih et 
alP 

The quadratic version of Eq.(2.1.3) together with the standard k-e transport 
equations, successfully predicts many complex flows as well as simple flows which 
include backward- facing step flows; confined coflowing jets; confined swirling coaxial 
jets; flows in 180° curved duct; flows in a diffuser and a nozzle; boundary layer 
flows with pressure gradient and turbulent free shear flows. See references 23-25 for 
detailed results. 

The Researchers involved with the subject in this section are J. Zhu and T.-H. 
Shih. 

3.4 Scalar flux algebraic equation model 

In parallel with Reynolds stress algebraic equation model, we have also tried to 
develop an effective anisotropic scalar eddy diffusivity model for scalar (heat) fluxes 
based on the new constitutive relation (2.1.7) and the new thermal dissipation rate 
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equation (2.2.16). We have determined that it seems sufficient to truncate Eq.(2.1.7) 
up to linear terms of the mean velocity gradient, i.e., 

Wi = + T(lT ]1/2{a * Ui ’ i+a3Ui ' i)T ’ i (3 ' 41) 

c tQ c 6 £$ 

This equation indicates that the heat flux and the mean temperature gradient are 
not necessarily in alignment due to the distortion of the flow field. This means that 
the effective scalar eddy diffusivity is anisotropic. 

Eq.(3.4.1) together with the 0 2 and eg equations will be a closed set of model 
equations for turbulent heat fluxes. The model coefficients are calibrated from 
homogeneous flows. Detailed analysis and a few model tests are described in this 
research briefs by A. Shabbir. 

The Researchers involved with the subject in this section are A. Shabbir and 
T.-H. Shih. 

3.5 Reynolds stress transport equation model 

The Reynolds stress transport equation model is considered as a next generation 
of advanced turbulence modeling for engineering applications. In principle, the 
second moment equations describe various effects of the mean flow and external 
agencies on the evolution of turbulence, hence, are the most attractive way (also 
the simplest correct way) to study turbulent flows. 

Various closure models for second moment equations have been developed. The 
success of these closures are marginal and vary with each flow. To identify the 
sources of their deficiencies, one often uses simple flows where the specific model 
term in the second moment equations can be isolated, hence, the corresponding 
model can be checked against experimental data or direct numerical simulation 
(DNS). For example, using pre-distorted anisotropic homogeneous relaxation flows, 
we may check the return-to-isotropy models with experimental data or DNS. How- 
ever, for other flows, several model terms, such as, triple velocity correlations, rapid 
and slow pressure- strain correlations, etc., simultaneously exist and can not be iso- 
lated in the experiments. In these cases (for example, in a homogeneous shear flow 
or a channel flow) only DNS can provide all the information for simultaneously 
checking various models. 

We have examined various existing closure models using experimental data as 
well as DNS data (Shih et al 26 and Shih and Lumley 27 ). Conclusive statements are 
difficult to draw at this time. However, the following remarks can be made about 
various closures for the second moment equations, i.e., the triple velocity correlation 
Tijk, the rapid and slow pressure related correlations n*J, nfj, and the dissipation 
rate tensor 

a) T^k • All the existing models, such as Daly and Harlow 28 , LRR 10 , Lumley 29 , 
etc., are not very satisfactory for highly inhomogeneous flows, such as flow near the 
wall. However, for flows where the inhomogeneity is not very high, the above closure 
models become close to each other and also closer to the DNS data. In addition, 
the triple velocity correlations in these situations are usually small comparing with 
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other terms in the equation, so that modeling of this term is not as critical as other 
terms for the results of turbulent flow calculations, except for the flow near the wall. 

b) ng\ It is very clear from all the available DNS data that nonlinear models, such 
as, Shih and Lumley 30 are much better than linear models, such as SSG 31 . It seems 
also that the following constitutive relation 

nj = FQmt'Uij) 

is quite appropriate, i.e., its dependence on turbulent Reynolds number and other 
parameters is quite weak and can be neglected. However, one deficiency of this form 
observed by Reynolds 32 is that it can not take the rotation effect into account. 

c) nfj. This term is usually modeled together with the dissipation tensor £ij and 
the combination of the two is called the return-to-isotropy term. All existing mod- 
els are unsatisfactory at the present time. They are far from “universal”, i.e., their 
performance varies from flow to flow. It is noticed that some strange behavior of 
return-to-isotropy (for example, for some pre-distorted flow relaxation, turbulence 
evolves toward anisotropy before it returns to isotropy) occurs and cannot be pos- 
sibly modeled with the following constitutive relation: 

Ilfj = F(uiUj,k,e) 

In addition, the behavior of return-to-isotropy was found to depend not only on the 
Reynolds stresses at the present time but also on their history according to DNS 
data (Lee 33 ). It may be also necessary to include triple velocity correlations into 
the above constitutive relations from the definition of n^. The term Ilf - seems 
highly dependent on the turbulent Reynolds number and slowly approaches to its 
asymptote as Reynolds number goes to infinity, so that, in general, one should not 
exclude its dependence on turbulent Reynolds number even for moderate Reynolds 
numbers. In addition, IIfj is also noticeably affected by the mean strain rate ac- 
cording to the DNS data 34 , so that, in general, the mean strain rate should be also 
considered in the constitutive relation. In short, much more research is needed for 
developing a better model of II*j. 

The Researchers involved with the subject in this section are T.-H. Shih and A. 
Shabbir. 

3.6 Non-equilibrium multiple-scale model 

To consider the effect of the nonequilibrium of energy spectrum on turbulent 
quantities, such as k, e and uiu], etc., Hanjelic et al . 35 are the first to propose a 
partition in the turbulent energy spectrum. Because of the nonequilibrium, the rate 
that energy enters the low wave number region, £ p , does not equal to the energy 
transfer rate from low wave numbers to high wave numbers, s t . Therefore it is 
reasonable to describe the evolution of the energy contained in low wave number 
region, k p , and high wave number region k t , separately. As a result, the time scale 
or the length scale defined by different energy transfer rates will be different and 
this multiple-scale concept reflects the nonequilibrium effect of turbulence. 
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We think that this concept would be more appropriate for compressible flows 
because the compressibility often creates nonequilibrium interactions between large 
and small eddies. We first modify Hanjelic et al.’s model, test it in various free 
shear flows and boundary layer flows and then extend it to compressible flows by 
consideration of the effects of compressibility on the equations for k p and e p . The 
proposed model is tested in both compressible free shear flows and boundary layer 
flows. For detailed analysis and flow calculations see the report by Duncan et al. Z6 
and Liou and Shih 37 . 

The Researchers involved with the subject in this section are W. Liou, T.-H. Shih 
and B. Duncan. 

3.7 Bypass transition model 

The onset of turbulence transition in the propulsion system is often highly influ- 
enced by the free stream turbulence. This transition process does not go through 
the linear instability but is mainly controlled by nonlinear processes. Therefore, it 
is sometimes called “bypass” transition. Because of this highly nonlinear process of 
transition, turbulence models may be used to predict it. In fact, many two-equation 
models, for example, k-e eddy viscosity models of Launder and Sharma 38 , Chien 15 , 
etc., do mimic bypass transition on a flat plate when the free stream has a certain 
amount of turbulent intensity. However, to obtain an accurate prediction of bypass 
transition, the study of the bypass transition process and physics is needed. The 
conventional turbulence models must be modified to take into account the intermit- 
tent phenomena of transitional flows. 

We have proposed transition models based on a two-equation turbulence model 
using an intermittency factor to modify either the eddy viscosity or modeled k-e 
equations. Successful results for a flat plate boundary layer under various free- 
stream turbulence intensities are obtained. For details see the report by Yang and 
Shih 39 . 

The Researchers involved with the subject in this section are Z. Yang and T.-H. 
Shih. 

3.8 Joint scalar PDF model 

One of the critical problems in turbulent combustion is how to treat the inter- 
action between the chemical reaction on the turbulence. The estimation of the 
production rate of compositions based on the mean flow temperature would be in a 
very large error for flows with finite rate chemical reactions. The reason is that the 
production rate of compositions depends not only on the mean values of tempera- 
ture T and compositions Ci, but also very much depends on the detailed fluctuations 
of temperature 0 and compositions c,. The moment closure scheme of modeling the 
production rate of compositions in terms of the mean flow temperature, the mean 
compositions and various correlations consisting of the fluctuating temperature and 
composition, such as 0 2 , Qa, cicj, ..., has not been successful. However, the PDF 
method allows us to treat chemical reaction exactly without modeling (Pope 40 ). 
Therefore, for the study of turbulent combustion problems, we use the joint scalar 
PDF transport equation for the scalar field and the moment closure schemes for 


Research Activities at CMOTT 


17 


the velocity field and develop a hybrid solver consisting of a Monte Carlo scheme 
and a conventional CFD method. For detailed description of this procedure and its 
applications see Hsu 41 and Hsu et a/. 42 

The Researchers involved with the subject in this section are A. Hsu, A.T. Norris 
and J.Y. Chent 

4. RNG and DIA 

In developing one point turbulence models, conventional modeling methods can 
be supplemented by “non- conventional” methods such as renormalization group 
theory (RNG) and the direct interaction approximation (DIA). These are two point 
theories formulated in wavevector or fourier space; one point models are derived 
by integration over wavevectors. This approach provides theoretical support for 
conventionally derived models and sometimes suggests theoretically derived forms 
for the empirical elements, whether constants or functions, which appear in these 
models. 

We have applied RNG methods to both the eddy viscosity and Reynolds stress 
transport equation models. In addition to the k — e model proposed by Yakhot 
and Orszag 43 , it is possible to obtain constitutive relations for Reynolds stress 
and heat fluxes (Rubinstein and Barton 4 ’ 6 which are special cases of the general 
results Eqs.(2.1.3) and (2.1.7). By applying the perturbation theories of Yakhot 
and Orszag 43 to the relevant correlations, expansions in powers of the mean velocity 
gradient are obtained for the stresses and heat fluxes; quadratic truncation of the 
series leads to a stress model Eq.(2.1.3) with constant a 4 ,a 6 ,a 7 and a heat flux 
model Eq.(2.1.7) with constant o 2 , a 3 in which the constants are in good agreement 
with empirically selected values. The forms derived are also consistent with the 
DIA analysis of Yoshizawa 3,5 . 

The RNG method also provides a formulation for closing the Reynolds stress 
transport equation (Rubinstein and Barton 44 ). Perturbative evaluation of the cor- 
relations n ^ and n?j leads to series expansions in powers of the mean velocity 
gradient. These series can be consolidated, or “resumed” using the known pertur- 
bation series for the Reynolds stresses by methods analogous to Pade approximation. 
Systematic lowest order summation leads to a Reynolds stress transport equation 
with a form identical to the LRR model equation and with constants in reason- 
able agreement with empirically chosen values. Higher order resummation leading 
to nonlinear models of the type described in Sec. 3.5 remains an open possibility. 
The possibility of such resummation in the context of DIA has been discussed by 
Yoshizawa 45,46 . 

Recent work has focussed on nonequilibrium time dependent relations between the 
Reynolds stress and the mean flow derived from a simplification of the DIA theory 
of shear turbulence. In this theory, shear turbulence is modeled by a non-Markovian 
eddy damping acting against the mean shear. The RNG and DIA Reynolds stress 
transport models and the LRR model all assume Markovian damping; as in the 
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molecular theory of transport coefficients, Markovian damping describes long time 
behavior and is incorrect at short times. The most important consequence of non- 
Markovian damping is a strong suppression of eddy damping at short times. This 
leads to closer agreement between the present theory and rapid distortion theory 
at short time, this is important in modeling oscillating shear flows: recent work of 
Mankbadi 47 shows that RDT based models best predict such flows. In transient 
homogeneous shear flow at high strain rates, the LRR model predicts rapid onset 
of eddy damping leading to excessive growth of turbulence kinetic energy at short 
times. The suppression of eddy damping at short times in the present model should 
lead to improved predictions for this flow as well. 

Another consequence of this theory is a stress model Eq.(2.1.3) in which the coef- 
ficients a 2 , . . . are functions of the mean strain rate. This theory can be described as 
RDT with a modified total strain determined by the response function of the DIA 
theory of isotropic turbulence. The introduction of a phenomenological modified to- 
tal strain has often been advocated in the RDT literature to improve the agreement 
between RDT and shear flow data; here the modified total strain is deduced as a 
consequence of the theory. In the special but important case of simple shear flow 
in which dUi/dxj = S 6 n 6 j 2 , the result can be formulated in terms of Eq.(2.1.3) in 
which, for example, 02 = a 2 (Sk/e ) and the function a$ is found exactly from RDT. 
There are analogous results for the coefficients a^,a^, 07 ] in simple shear flow, the 
remaining terms in Eq.(2.1.3) identically vanish. Extension of this theory to other 
mean shear tensors depends on the tabulation of the corresponding RDT solution. 

The researchers involved with the subject in this section are R. Rubinstein and 
A. Yoshizawa.f 

5. Numerical Simulation 

To obtain a better understanding of the effect of compressibility and rotation 
on turbulence, numerical simulations of compressible homogeneous shear flows and 
rotational flows are carried out. The effects of compressibility and rotation on the 
energy spectrum and energy cascade between turbulent eddies has been analyzed 
(Hsu and Shih 48 ). These simulations support the idea of the multiple-scale model 
for nonequilibrium compressible turbulent flows (W. Liou and Shih 37 ). 

Another numerical simulation is the transition subjected to the free stream large 
disturbances. The objective of this simulation is to obtain some insight into the 
transition physics and to provide data base for bypass transition modeling. Based 
on the assumption that the transition process is mainly controlled by large scale 
motions, we use a high accuracy finite difference Navior-Stokes solver with course 
grids to simulate the large scale motions of transition. A preliminary calculation of 
bypass transition was carried out. Various statistics of the calculated flow field are 
under examination. 

The Researchers involved with the subject in this section are A. Hsu, C. LiouJ, 
Z. Yang, A. Shabbir, T.-H. Shih. 

f Professor, Tokyo University, Japan 
| Professor, University of Colorado, Denver, CO 
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Research on Two Equation Models 
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Motiv ations and Objectives; 

• k-e model is the most widely used turbulence model in engineering calculations. 

• However, the following deficiencies need to be fixed: 

- Currently, most k-e models for near wall turbulence contain geometry parameters. 

— The form of the e equation lacks a solid theoretical foundation. 

-The k-e model performs rather inadequately for flows with adverse pressure gradient. 

- The capability of k-e models in predicting bypass transition due to the freestream 
turbulence needs improving. 


Modeling of Near Wall Turbulence 


• Near wall k — e model = Standard k — e model + near wall effect. 

• The near wall effect: 

— The time scale approaches the Kolmogorov time scale near the wall. 

- The damping function is parametrized by a new parameter which is independent of the 
coordinate system. 

• The resulting model is Galilean and tensorial invariant. 

• The resulting model is robust numerically. 













On the Wall Functions 


• The advantages of the wall function approach: 

— Reduce the number of grid points by half, at least. 

- Reduce the numerical stiffness of the dissipation rate equation, by less grid stretching. 

• The limit of wall function approach: the flow is assumed to be attached to the wall. 

• Existing wall functions axe based on the flat plate BL at zero pressure gradient. 

It is inadequate when the pressure gradient is not zero. 

• A new set of wall functions are obtained: 

- They are based on the asymptotic behavior of the governing equations in the log layer. 

— They contain the effect of the pressure gradient. 
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A Vorticity Dynamics Based Model for the Dissipatio n Rate Equation 


• The dyn ami c equation for the fluctuating vorticity is analyzed. 

(The terms in the fluctuating vorticity equation have clearer physical meanings than terms 
in the dissipation rate equation.) 

• For large Reynolds numbers, e = vts^o\. 

• The resulting model equation has a better foundation than the standard e equation. 

• The resulting model equation always gives, a positive production in dissipation rate. 

- The model calculation is expected to be more robust for complex flow calculations. 


Flow Inhomogeneity and the Dissipation Rate Equation 

• The exact dissipation rate equation contains source terms due to the flow inhomogeneity. 

• However, the existing e equations are homogeneous. 

(The source terms are the same for both the homogeneous and inhomogeneous flows.) 

• A new model equation for e is proposed, which accounts for the inhomogeneity effect. 

— Flow inhomogeneity is represented by V5 and Vfc. 

— Invariant theory is used to derive such a model equation. 

• The resulting model equation accurately account for the effect of the pressure gradient. 
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Mode ling of Bypass Transition 


• Low Reynolds number k — e models could mimic transition. 

- The predictions are not very good. 

— Among these models, the Launder-Sharma model gives the best prediction. 

• New model for bypass transition is proposed by introducing the effect of the intermittency 

• The transition model recovers to turbulence model at the end of the transition zone. 

• Calculations of the benchmark flows show that the present model gives better predictions 
compared with the k — e models without the intermittency effect. 
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.S ummaries: 


The capabilities oik - e model are enhanced in the following areas: 

• A Galilean and tensorial invariant k — e model for near wall turbulence. 

• A new set of wall functions for attached flows. 

• A new model equation for the dissipation rate: 

- It has a better theoretical basis. 

- It contains the contribution of flow inhomogeneity. 

- It captures the effect of the pressure gradient accurately. 

• A better model for bypass transition due to freestream turbulence. 



A REYNOLDS STRESS 
ALGEBRAIC EQUATION MODEL 


Presented by 
J. Zhu 


Motivation 

• Two-equation turbulence models (K-e, K-o;) 

- Simple, robust and 
computationally inexpensive 

- Isotropic 

- Not realizable, u jui may become negative 

• Reynolds stress transport equation models 

- Anisotropic 

- Capable of simulating u^ij-transport 

- Difficult to model higher-order correlations 

- Complex and computationally expensive 



Objective 


Develop a Reynolds stress algebraic equation model 
which combines the respective advantages of 
two-equation models and 
Reynolds stress transport equation models: 

• Simple, robust and efficient 

• Anisotropic 

• Realizable 


Derivation of Constitutive Relation 

• First step: assume some relationship between 
quantities of interest v^uj = Fij(Uij, K, e) 

• Second step: use invariant theory to determine 
the function F{j 

1. An invariant can only be a function of other 
invariants 

2. If the LHS of an equation is bilinear in arbi- 
trary vectors A{ and Bj, the RHS must also 
be bilinear in and Bj 


• A general constitutive model of uiUj 


u Wj = |*«B + 2e^(Uij + Uj, , - | Ui.Ai) + 2 <n§r(°y + ^ 

+ 2a^~-(Ui l kUj l k - j I + 207 -^(Uk,iUkj - jn 2 5y) 

+ la^V^kUlk + UlkVj* - |n,« B ) + 2 « 1 o^(Di, i O?j + «X. - | n 3«„) 
+ 20,^0 ilkVlk - in.%) + - 5 na{ ‘i) 

+ 2<*u^-(Ui,kV,, k u?j + Uj.kU^kUli - ^n 5 <,j) 

7^6 2 

+ 2ai 6 — (Ui t kUl k Ulj + UjjUfcUt'i - -n 6 5jj) 

+ 2ax8 (Ui^Ui^k U^ m Uj, m + Uj^kU^U^ — -IMy ) 


n x — Ui, k u kyi , 
n 5 = u^u^ul , 


n 2 = Ui t kUi ik , n 3 = , n 4 = u£ k uf tk , 

lie = u itk u? tk u? ti , n 7 = Ui t kU[ t kUi tTn u i rn 


The first five terms are of the same form as that of Yoshizawa (1984) 
and Rubinstein and Barton (1990). 


• Realizable Algebraic Equation Model 


u i u j = ^ K6ij - v t (Uij + Uj } i) 


. C Tl K 3 
A 2 + V 3 + £ 3 e 2 


(U ijk U k j + U jjk U k j - 



Ct2 if 3 

A 2 + tj 3 + £ 3 e 2 

C T 3 if 3 

A 2 + *? 3 + £ 3 e 2 


{Ui, k Uj,k - \ftkj) 
{UkjUkj - 


K t t + [UjK — (y + — u i u jUiJ e 


«,t + Wjt - O + — ~ c lj( u i u j u i,j c 2 k 



where 


2/3 #2 
Vt Ai + t) + € 

>? = ^, S = (2S ij S l j) 1 / 2 , S;,- = |(£/,,J + Uj,i) 

( = ^p-, n = ( 2 n; j n‘ } ) 1 / 2 , = |<cr iJ -^ j )+4 £mijUro 

C Tl = -4, C r2 = 13, C r3 = -2, Ai = 5.5, a = 0.1, A 2 = 1000. 
Ci = 1.44, C 2 = 1.92, a K = 1, a e = 1.3 


Solution Procedure 

• Conservative finite-volume method 

• Non-staggered grids 

• Momentum interpolation for velocity-pressure 
coupling 

• Second-order accurate differencing schemes 

• Quadratic terms in v^u] treated as source 

• Strongly implicit solution algorithm of Stone 


Performance in Complex Turbulent Flows 


Diffuser flows 

Backward-facing step flows 
Confined jets 

Confined swirling coaxial jets 


Diffuser flows 

• Fraser’s case, a = 10° 

• Trupp, Azad & Kassab's case, a = 8° 


Wall pressure 



SKE : Standard K-e model 
RNG : RNG-based K-e model 
RRSAE : Present model 
KO : K-uj model 
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Backward-facing step flows 


• Driver Si Seegmiller’s case, H s /H d = 1/8 

• Kim, Kline Si Johnston's case, H s /H d = 1/2 



Comparison of reattachment points 
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Reynolds stresses 
(DS-case) 



Confined jets 

• Barchilon & Curtet's case, 6 — 0 ° 

• Binder & Kian’s case, 6 = 2.5° 































Confined swirling coaxial jets 


• Roback & Johnson’s case 
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Conclusions 


• A general constitutive relation for UiUj has been 
derived 

1. its quadratic form seems sufficient 

2. It is realizable 

• It significantly improves the performance of K-e 
based models 

• Easy to implement 

• No significant increase in computing time 


• Very robust 



ASSESSMENT AND DEVELOPMENT OF SECOND ORDER 


TURBULENCE MODELS 
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Presented by 
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MOTIVATION 

• these models describe the effect of mean flow and external 
agencies (such as buoyancy) on the evolution of turbulence 

• therefore, in principle, these models give a more accurate de- 
scription of complicated flow fields than the two equation mod- 
els 

- e.g flows with large anisotropy in turbulence (such as near 
the leading edge of a turbine blade) 
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OBJECTIVE 


• assess the performance of the various second order turbulence 
models in benchmark flows 

• seek improvements where necessary 

- model for the pressure correlation term in the jscalar flux 
equation 

- model for the scalar dissipation equation 



Transport Equations for Second Moments 


DuiUj 

ij 


+ Uij + Tij — D{j 


Du t 6 

Dt 

d¥ 

~Dt 


=Pi + IT + Ti 
=P + T — D 


These equations have to be closed by providing models for: 

• Pressure correlation terms (Hj, 11*) 

• Transport (Diffusion) terms Ti) 

• Dissipation terms (Aj, D) 


3 


HOW TO ASSESS MODELS ? 


Global computation 

• Mean and turbulence equations are numerically solved 


Du i u j =P .. + 
Dt 13 


• Results (e.g. Reynolds stresses) are then compared with ex- 
periments or DNS data 
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Direct comparison 

Individual terms in the turbulence equations (such as pressure cor- 
relation terms) are directly compared with experiment or DNS data 

Note that: 

- In experiments pressure correlation terms can not be mea- 
sured but can only be obtained indirectly through balance 
of second moment equations 

- DNS allows direct computation of these correlations but 
is limited to low Reynolds number 

Most of the results to be shown in this presentation are direct com- 
parisons 
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Models for [pressure correlation term in the scalar flux equation 


• a. 

• b. 


• c. 

• d. 


Launder model (1975) 

Zeman and Lumley model (1976) 

- linear in scalar flux 

- do not satisfy realizability 

Shih, Shabbir and Lumley model (1985,1991) 
Craft, Fu, Launder, Tselepidakis model (1989) 

- linear in scalar flux and Reynolds stress 

- satisfy realizability 
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Application to Homogeneous Shear Flow 
(Experiment as well as DNS data) 
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CONCLUSION 


Models for pressure correlation term in scalar flux equation 

• Models involving both scalar flux and Reynolds stress give bet- 
ter performance than the models which involve only scalar flux. 
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Models for pressure correlation term in the Reynolds stress equation 

• a. Launder, Reece, Rodi model (1975) 

• b. Speziale, Sarkar and Gatski model (1991) 

- linear or quasi-linear in Reynolds stress 

- do not satisfy realizability 

• c. Shih and Lumley model (1985) 

• d. Craft, Fu, Launder, Tselepidakis model (1989) 

- nonlinear in Reynolds stress 

- satisfy realizability 
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Tavoularis & Corrsin Expt (1981) Tavoularis & Corrsin Expt (1981) 













DNS Data of Rogers et al. (1986) DNS Data of Rogers et al. (1986) 






















CONCLUSION 


Models for pressure correlation term in Reynolds stress equation 

• For the DNS data non-linear models give better performance 
than linear models. However, for the experiment no single 
model performs better for all the components 

• For the rapid part of the pressure correlation the relation 1^ = 
F(uiuj, U{j) is found to be adequate 
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CONCLUSION (contdj 

• Performance of all the slow pressure correlatiion models varies 
from one flow to another. 

• Furthermore, the relation Lt^ = F(uiu],k , e), is inadequuate 
in certain situations 

- DNS shows that II - is dependent not only on the present 
time value of Reynolds stress but also on its past history 

- Definition of 11^- implies that it is also a funciton of triple 
velocity moment uiujuk 

• Therefore, more research is needed before any model for 1^- 
can be recommended for use. 
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A New Model Equation for Scalar Dissipation 

• Traditional scalar dissipation rate equation is modeled in an 
analogue fashion to the mechanical dissipation equation 

• Equation proposed here is modeled after the exact equation 
for scalar dissipation 

• Its production/destruction mechanisms are different than the 
traditional model quation 
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Application to Homogeneous Benchmark Flows 

1. Homogeneous turbulence subjected to constant scalar gradient 

2. Homogeneous turbulence subjected to constant scalar gradient 
and constant shear 

Global computation of the following two equations 


d6 2 — -dT n 

Ui dT j ^~ ' 9 dl % 6 

V °JL 

’dx. 


— Coi€oS + Cq2 


\fPr 


-C e 3 


606 

~k 


Mechanical field (i.e. fc, e, etc.) and scalar flux,iii#, are taken as 
known. This way performance of the scalar dissipation equation is 

isolated. 
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Tavoularis & Corrsin Expt (1981) 












CONCLUSION 


• The transport equation for thermal dissipation rate proposed 
here gives improvement over the standard equation in at-least 
all the simpler benchmark flows. Its performance in the wall 
bounded flows is being assessed 
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A Multiple-Scale Model for 
Compressible Turbulent Flows 


presented by 


William W. Liou 


Workshop on Computational September 15-16, 1993 

Turbulence Modeling 


Outline 

□ Motivation 

□ Model Development 

• formulation and determination of model coefficients 

□ Model Predictions: 

• compressible shear layer 

• compressible boundary layer 

□ Concluding Remarks 


Motivation 


□ Incorporate the effects of non-equilibrium energy spectrum 
of compressible turbulence in computational models 


DNS of Compressible Turbulence 

□ formation of high gradient regions or eddy shocklets 

□ enhanced vortex stretching and spectral energy transfer 

□ non-equilibrium turbulent kinetic energy spectrum due to 
the flow compressibility 


Multiple-Scale Model 



Model Equations 

□ Large Scale 



d r/ _ 

+ <7 


+ M d i? 


dy' 


p€p + fCi 



3 r/— - \ dtp -| , 

= ^ [(f,+ ^F ) ^ 1+ 


CpiS-Ar (^) 2 

°y 


~ 2 

Cp2p~^~ + f c 2 

k v 


• fc 1 — exchanges between the turbulent kinetic energy and internal energy 

• fc 2 - increased spectral energy transfer due to eddy-shocklets 


Model Equations (cont.) 
□ Small Scale 


_Dk t 

? dT 


_De~ 
9 Dt 


= > + ^f 1 + CW f- CW f 


□ Eddy Viscosity 


HT P U l ~ p[kp + fci} 5 


i (kp + kt)+ 


Model Coefficients 

□ Decaying turbulence and homogeneous turbulence 


Cpi = (1 - 2) + ^Cj>2, Cp 2 
a a 

Cti = 1 - i Ct 2 


71+1 


n 


(3-1 + C P2 (3f 

Kp 

(3 + (3&- 1 

kp 


C, 

% 

(Jr 

k t 



n 

a 

/? 

0.09 

d5 — 
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Model Equations (cont.) 

□ Large Scale 

• fci - Sarkar et al. (1992) 

&VL 

f Cl = -a:2Mt/Lir(^) 2 + «3 M?pe p 

• fc 2 - dimensional analysis 

fC2 — f c 2 (p? Pj € > p) 

p i n = ]*L = — 

ni = ~ tjwt 2 ^ 

fc 2 = 6i + 6 2 M, 2 + 0(M t 4 ) 

Model Predictions 

□ Compressible Free Shear Layers 



• mean velocity profiles. M c =0.51 

• peak Reynolds shear stresses. M c =0~l.6 

• vorticity thickness growth rates 









Multiple-Scale Model 
□ Wall Function 




Model Predictions 


□ Compressible Boundary Layers 



• mean velocity profiles. Moo=2.831, Re$=420,100 

• mean temperature profiles. 

• skin friction coefficients. Moo=0~5, it*e 0 =lO 4 










Concluding Remarks 

□ The present compressible multiple-scale model predicts 
correctly both the spreading rate of compressible shear 
layers and the skin friction coefficient of compressible 
boundary layers. 

□ Need to implement the model into the calculation of more 
complex compressible turbulent flows. 
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CLOSURE PROBLEM: 


Ui = Ui + u\, 

¥> = ¥> + y/, 


u'-u'j — Turbulence Modeling 

Y"u" — Analogy of shear stress: diffusion model. 
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ANALYTICAL THEORIES OF TURBULENCE 
APPLIED TO TURBULENCE MODELING 



presentation by 
Robert Rubinstein 
CMOTT 


MOTIVATION 


Turbulence models contain undetermined elements 

constants 
functions 
model forms 


Apply analytical theories to determine these elements 

direct interaction approximation (DIA) (Kraichnan, 1959) 
renormalization group (RNG) (Yakhot and Orszag 1986) 


EXAMPLE 1 

Determine constants in nonlinear eddy viscosity (NLEV) (Yoshizawa 1984) 
and nonlinear eddy diffusivity (Yoshizawa 1987) models 


< uu > -hi = C„— (VU + VU T ) 

o ^ 

+ C t1 ^(VUVU t - Ivu : VU T I) 

+ C T 2-y(VUVU + VU T VU T - ^VU : VU I) 

e 2 £ 

+ CV 3 ^(VU T VU - ivU T : fU I) 

< <?u > = cXve + o^vu • w + c« 2 ^v u T • ve 


Universal inertial range constants for isotropic turbulence. 

E(k) = C K e 2/3 k- 5/3 
u(k) = Cl)€ 1/3 fc- 4 / 3 

C K = 1.61 C D = 0.49 (Yakhot and Orszag 1986) 

Assume shear is a weak perturbation of isotropic 
background state. Then (Rubinstein and Barton 1990): 


exact 

value 

experiment 

C, (\?CdI<% 

0.085 

~ 0.09 

Cr 1 lfe(|) 2 /(Cl>C K ) 2 

0.034 

~ 0.04 

c r3 nK§) 2 /(cw *) 2 

-.014 

~ — .01 


Passive scalar results 

Pr? = f(Pr, Re) (Yakhot and Orszag 1986) 

/ is a known function: for Re — oo, Pr? ~ 0.7 
C K = CJPr T 

Constants of the nonlinear theory are known fimctions of Ct>, CkiP t t 
(Rubinstein and Barton 1991) 

Diffusivity ratio * 42/^22 = 1.9 (Tavoularis and Corrsin), 2.3 (theory) 

In principle, higher order nonlinearities in NLEV 
can be computed theoretically 


EXAMPLE 2 

Determine model form for rapid pressure strain correlation II 

RNG gives a series for II in powers of VC/. 

The LRR model 

II = Ct\ < uu > VU + C+2 < uu > VU T 

arises from systematic lowest order consolidation (resummation) of this 
series. (Rubinstein and Barton 1992): 



LRR 


RNG 

Cr 1 

0.7636 

16 

21 

- 0.7619 

CV2 

0.1091 

2 

21 

~ 0.0952 



Second order resummations are possible: 
II = F(< uu >)VU + G(< uu >)VU T 


n =< uu > [Ci vu + c 2 vuvu + • • ■] 


These are speculative at this point. 


EXAMPLE 3 

Determine unknown functions in generalized nonlinear eddy vis- 
cosity representation for simple shear flow 

Nonlinear eddy viscosity does not apply at large strain rates: it predicts 
< uu > /k — ► oo. Introduce the generalization 

< uu > -hi = CM — (vu + vu T ) + • • • 

o £ 

where in simple shear flow dUi/dxj = SSnSj2, V — Skfe 


C-3 



DIA analysis of shear turbulence leads to 


^ U{Uj (^Tj j G R^j 


where F l3 is known from rapid distortion theory (RDT) (Maxey 1982) and 
Cr is a universal constant ( Cr ~ 1.59, Yakhot and Orszag 1986). 

General time dependent theory 


< U{Uj > ( t ) = kFij(a*(t )) 


where a* is a modified total strain determined by DIA: 


«•(*) 



ds G(t — s)'S(s) 


G is the response function of DIA. 


Homogeneous shear flow at high strain rate 
T) = Sk/e — * oo 



& 12 

bn 

b 22 

b 33 

NLEV 

0(7?) 

0( v 2 ) 

o(v 2 ) 

0( T? 2 ) 

LRR 

0( v- 1 ) 

f: 11 ( c ) 

/22(c) 

/33(c) 

present 

0(77 1 ) 

2/3 

-1/3 

-1/3 

RDT 

0(t7 -1 ) 

2/3 

-1/3 

-1/3 


In LRR, fij(c) are model dependent constants. The limiting normal stresses 
are not necessarily realizable. 

DNS (Lee, Kim and Moin) shows that RDT is a good description of these 
flows even for moderate times. 



CURRENT WORK 


Time dependent model for nonequilibrium turbulence based on DIA 

New evaluation of NLEV constants 

Nonmarkovian damping: short time suppression of eddy damping. 

Required for oscillating and highly strained flows. 

Two scale theory of shear turbulence: 
large scales RDT, 

small scales ~ Kolmogorov inertial range dynamics 
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DNS for Bypass Transition: Cases Calculated 


• 2D/3D linear instability stage (temporal simulation) 

• 2D/3D linear instability stage (spatial simulation) 

• Secondary instability stage (spatial simulation) 

• Simulation of the whole transition zone (spatial simulation) 



comparison of streamwise distu rbance velocity 
DNS+MC 





DNS: Effect of Rotation 



EQUATIONS FOR TRIPLE CORELATION EQUATIONS 
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DNS: Reactive Flows 







£4j*tfeddy*= 7.44E-01 qO= 5£4E«0I 
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Title 


ALGEBRAIC MODELS FOR TURBINE 
BLADE HEAT TRANSFER 

R. J. Boyle 

NASA Lewis Research Center 


Workshop on Computational Turbulence Modeling 


LeRC j 9/93 Aim and Scope [_2 

• VALIDATE PREDICTIONS OF TURBINE 
BLADE HEAT TRANSFER & PRESSURE LOSS 

SMOOTH AND ROUGH BLADES 

• ALGEBRAIC - BALDWIN-LOMAX TYPE 
MODELS 

• IMPLEMENTED IN NAVIER-STOKES CODES - 
RVC3D & RVCQ3D 

• ENGINEERING APPLICATIONS 


Workshop on Computational Turbulence Modeling 











SMOOTH SURFACE PREDICTION 
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• FREESTREAM TURBULENCE - ex LEADING 


EDGE 

• TRANSITION MODEL - ex MAYLE'S 

• ALGEBRAIC - BALDWIN-LOMAX TYPE 
MODELS 


Workshop on Computational Turbulence Modeling 
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TUoo EFFECT ON LEADING EDGE HEAT TRANSFER 



Workshop on Computational Turbulence Modeling 






Workshop on Computational Turbulence Modeling 
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APPROACHES FOR ROUGH SURFACES 

• TAYLOR. COLEMAN, and HODGE 

EXPLICITLY ACCOUNT FOR BLOCKAGE, 
DRAG AND HEAT TRANSFER 

SOURCE TERMS IN NAVIER-STOKES 
EQUATIONS 


• CEBECI and CHANG - USED IN ANALYSIS 
MODIFY MIXING LENGTH 
ONLY EDDY VISCOSITY AFFECTED 


Workshop on Computational Turbulence Modeling 
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Workshop on Computational Turbulence Modeling 
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CORRELATION FOR EQUIVALENT HEIGHT RATIO 
FROM SIGAL and DANBERG 


EQUIVALENT 
HEIGHT 
RATIO, h eq /h 


\ \\ 

\ 

\ \ \ 


\ \ \ 
V v ' 


\ \ \ \ 

\ \ \ \ 

\ \ \ \2D 

\ \ ^3D SIGAL & DANBERG 

\ \ 

\ \ 

DIRLING\ \ DVORAK 


100.0 


1000.0 


ROUGHNESS DENSITY PARAMETER, A 


Workshop on Computational Turbulence Modeling 
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EXPERIMENTAL DATA USED FOR COMPARISONS 

• HEAT TRANSFER 

HOSNI - FLAT PLATE 

BLAIR & ANDERSON - LARGE LOW SPEED 
ROTOR 

DUNN & KIM - SSME FUEL TURBINE 

• TURBINE EFFICIENCY 

BOYNTON - TWO-STAGE SSME HIGH PRESSURE 
FUEL TURBINE 

Workshop on Computational Turbulence Modeling 



Workshop on Computational Turbulence Modeling 
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RATIO OF PREDICTED TO MEASURED HEAT TRANSFER 

FOR A FLAT PLATE. Re x =400,000. 

(DATA OF HOSNI et al.) 



Equivalent Height Correlation 


Sigal &Danberg 


hoq/h Stp/St ro 



Dvorak 

^eq/h 

Stp/Stjn 

2.43 

0.50 

0.061 

1.85 

1.57 

1.15 


Dirling 


1.69 

0.45 

0.079 



1.10 


Workshop on Computational Turbulence Modeling 
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PREDICTED AND MEASURED BLADE HEAT TRANSFER 
(Data of Blair and Anderson) 

10 1” □ ROUGH A 

O - SMOOTH Ppl 


STANTON 

NUMBER 

X 1000 


SUCTION 


PRESSURE 


uU \ □ □ 

0 - 0 - 0 '°''° 


n ^ 1 t 

D o l ,r 

O Q 1 • 

' O O O ' 1 X 

o' o ^ u v ; 1 m □ 

\ 1 


1 c — -r n 

^ooocroo 


- 1.0 - 0.5 0.0 0.5 

SURFACE DISTANCE, s/b x 


Workshop on Computational Turbulence Modeling 

















Workshop on Computational Turbulence Modeling 




Workshop on Computational Turbulence Modeling 
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ERROR DUE TO SURFACE RESISTANCE 



TEMPERATURE MEASUREMENTS BENEATH ROUGHNESS 
AFFECT HEAT TRANSFER COEFFICIENT 


^ACTUAL / ^MEASURED = 1 + 

h - HEAT TRANSFER COEFFICIENT 

t - ROUGHNESS THICKNESS 

k - THERMAL CONDUCTIVITY OF ROUGHNESS 

• POSSIBLE REASON FOR HIGHEST AUGMENTATION 
CORRESPONDING TO LARGEST SCALE 


Workshop on Computational Turbulence Modeling 
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REQUIREMENTS FOR BENCHMARK ROUGHNESS 
HEAT TRANSFER DATA SETS 

• COMPLETE DESCRIPTION OF THE ROUGHNESS 
CHARACTERISTICS 

• SURFACE NOT DISTURBED BY MEASUREMENT 
TECHNIQUE 

• KNOWN TEMPERATURE OF ROUGHNESS ELEMENTS 

• VERIFICATION OF ACCURACY BY MEASUREMENTS ON 
SMOOTH SURFACES 


Workshop on Computational Turbulence Modeling 
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• SMOOTH SURFACES 

REASONABLE ACCURACY 
TRANSITION BEHAVIOR GOVERNS 
ACCURACY 

• ROUGH SURFACES 

INCONSISTENT ASSUMPTIONS 
NEEDED FOR LOSS AND HEAT 
TRANSFER 


Workshop on Computational Turbulence Modeling 
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NAVIER-STOKES TURBINE HEAT TRANSFER 
PREDICTIONS USING 


TWO-EQUATION TVRBULENCE CLOSURES 
ALI A. AMERI 

CENTER FOR RESEARCH, INC. 
UNIVERSITY OF KANSAS 
NASA LEWIS RESEARCH CENTER 


Internal Fluid Mechanics Division 
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TURBINE HEAT TRANSFER 


OUTLINE 


• Brief description of the method of solution 

• 2-d examples using various models 

• 3-d calculations using algebraic model and q-(£> model 


• Conclusions 




TURBINE HEATTRANSFER 


TURBULENCE MODELS 


LOWREYNOLDS NUMBER 

• COAKLETS c/-co MODEL 

• CHIEN’S K-z MODEL 


ALGEBRAIC MODEL 

• BALDWIN -LOMAX - WITHOUT TRANSITION MODEL 



Internal Fluid Mechanics Division 



TURBINE HEAT TRANSFER 


FORMULATION 


. MASS AVERAGED COMPRESSIBLE N-S EQUATIONS 
. + MODEL EQUATIONS 
• NON-PERIODIC C GRID. 




SOLUTION OF THE NA VI ER -STOKES EQUATIONS 


• FINITE VOLUME DISCRETIZATION 

• 4 STAGE RUNGE-KUJTA SCHEME 

• EIGENVALUE SCALING OF ARTIFICIAL DISSIPATION 

• VARIABLE COEFFICIENT RESIDUAL SMOOTHING 

• MULTIGRIDING 


Internal Fluid Mechanics Division 


TURBINE HEAT TRANSFER 


SOLUTION OF 'THE MODEL EQUATIONS 


• FINITE VOLUME DISCRETIZATION 

• 4 STAGE RUNGE-KUTTA SCHEME 

• VARIABLE COEFFICIENT RESIDUAL SMOOTHING 


• MULTIGRIDING 
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2- D CASES 

• SSME HIGH-PRESSURE FUEL TURBINE 1ST STAGE VANE 

• ALLISON’S C3X VANE DATA OF HYLTON ET AL. 

3- D CASES 

• LANGSTON’S CASCADE, THICK & THIN INLET BOUNDARY 

LAYERS. 


Internal Fluid Mechanics Division 


TURBINE HEAT TRANSFER 



Surface Pressure Distribution- First Stage Stator 









TURBINE HEAT TRANSFER 



s/c 


First Stage Vane, High Reynolds Number 



TURBINE HEAT TRANSFER 


SSME HP. FUEL TURBINE 1ST STAGE VANE 
Re 2=3. 6E5 


□ DATA 

q-w Mode! 

Chien Mode! 

BAL-LOM 

Laminar 


Pressure Surface 


Suction Surface 


First Stage Vane, Low Reynolds Number 


-T -- _J «. .. i.. . 1 _ r _> • . . - - .1. . ■ ■ Jt .• j.- a.-. } •=.. . ■. . •... ii - J . . i.- ' . -• 

1 - V' i 

Internal Fluid Mechanics Division 




TURBINE HEAT TRANSFER 


: Unrit 'Cftiter • 


SSME HP. FUEL TURBINE, 1ST STAGE VANE 
Re_2=3.6E5 

VARIATION OF LENGTH SCALE, q-w MODEL 


\ D 


AT/ 
1=0.010x0 
1=0.025x0 
1=0.050x0 


Pressure Surface 


Suction Surface 


First Stator, Variation with length scale 
q - co Model 








S/C 


First Stator, Variation with length scale 
Chien s Model 


C3x 1 49 

392x64 calculation, dy=1.e-5 



- 1.00 


cutsnr. 


10/04/91 13:17:22 




Thick Inlet Boundary Condition 


0 *- 



Experiment Graziani et.al. 




Computation, TRAF3D 



Langston’s Cascade 
Endwall Stanton Number xIOOO. 



Experiment, Graziani et al. Computation, TRAF3D 

Thick inlet boundary layer. 


THICK BOUNDARY LA YER INLET CONDITION 



q-omega Model, 193x49x65 


Experiment Graziani 





Surface Stanton Number x 10 00. , thin inlet boundaryjayer 






Langston’s Cascade 
Endwall Stanton Number xIOOO. 



Experiment, Graziani et al. Computation, TRAF3D 

Thin inlet boundary layer. 


LANCS ION’S CAS< 'ADH 
ENDWALL STANTON Nl iMliERx 1000 



Experiment, Graziani cl. al. M (n hkhIcI 193x49x65 grid 




TURBINE HEAT TRANSFER 


CONCLUSIONS 


• Two-equation models perform well in the fully turbulent regime. 

• The uncertainty re. length scale ofturb. needs to be resolved. 

• Transition modeling is crucial to the success of H-T analysis. 

• Need further improvements in convergence speed. 
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Objective and Rationale 


Objective: Improved prediction of turbine blade surface 
temperature and heat-transfer distributions 


Rationale: Existing heat-flux models lack the fundamental 
generality required for complex flows 




LeRC 


9/93 


Approach 


3 


Fundamental analysis of thermal DNS results 


Detailed evaluation of existing heat-flux models 


Development of thermal time scale model 
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Kasagi, Tomita, and Kuroda DNS 


5 



Pr . / 0.9 T e / T i 0.5 

turb 

Pr 4 constant T e / T ^ constant 

turb 


LeRC 9/93 Evaluation of existing heat-flux models (s> 


Constant Pr^ models inadequate 


Flux-transport models expensive 


Two-equation thermal models offer affordable improvement 
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Two-equation thermal models 


Hattori, Nagano, and Tagawa (HNT) 

complex damping functions (evaluation in progress) 


Sommer, So, and Lai (SSL) 

Shima coincidence condition (evaluation in progress) 


Thermal time scale model (development in progess) 
simple wail boundary condition: x Q = o 
simple near-wall asymptote: x ~ y 2 / 2 a 
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A COUPLED IMPLICIT SOLUTION METHOD FOR 
TURBULENT SPRAY COMBUSTION IN PROPULSION SYSTEMS 


K.-H. CHEN 

The University of Toledo, Ohio 
and 


J.-S. SHUEN 

Sverdrup Technology, Inc., Ohio 


Aerothermochemistry Branch 
Internal Fluid Mechanics Division 
NASA Lewis Research Center, Cleveland, Ohio 


OBJECTIVES 


• Develop an efficient and robust algorithm for multi-phase 
chemically reacting flows at all speeds, with emphasis on 
low Mach number flows. 

• Calculate turbulent spray combustion flow in a gas tur- 
bine combustor. 



MOTIVATION 


Many reacting flows in propulsion devices cannot be efficiently calculated by 
modern compressible flow CFD algorithms, e.g., 

rocket motor — wide range of Mach numbers, from near zero 
velocity at closed end to supersonic at nozzle exit, 
gas turbine combustor — low subsonic velocity, but large density 
variation precludes incompressible approach. 

Most low-speed reacting flow codes based on TEACH-type technologies - 
inefficient and lack of robustness for complex flows. 


Tremendous progress made in high-speed compressible flow CFD in past 
two decades. Extending application range to low-speed flow regime highly 
desirable. 


OUTLINES 


• GOVERNING EQUATIONS 

- Gas-Phase Equations 
— Liquid-Phase Equations 

• NUMERICAL algorithm 

• NUMERICAL TEST RESULTS 

• CONCLUSION 


FUTURE PLAN FOR ALLSPD CODE 


GOVERNING EQUATIONS 


• Gas-Phase Equations 


dQ dQ d(E - E.) 0(F - F.) 


r — + — + 

dr * + dr ^ 


r 

Q = — 
x J 



1/0 0 0 0 0 0 

u/P p 0 0 0 0 

v/P 0 p 0 0 0 

h/P — IpupvpOO 
k/P 0 0 ■ p • 

e/p 0 0 ■ • p 

Yi/P 0 0 0 ■ • 

Y 2 /p 0 0 0 0- 


Y N -i/P 0 0 


0 

2 c 3(^«) 
" 3 ° dx . 


S\p - Tge - 

+ a 3r 1 l 

y^-pe)-lSp t v(t + T y ) 
y s [(ci/l^ — C2fipli)^ + A] — 1 6ptv(fa + 3^) c i/1k 
y s S i 


y s 


E P n P^P ^ 

Ep UpThpUp 3" Pp r p n p it 

Ep flpfhpVp 3" Pp r p n p it 

Ep n p m v hj, — 4xr v n p hAT 
0 
0 

Ep«p>hp 


0 



ALLSPD MATN FEATURES 


Present Capabilities 

• 2-D and Axisymmetric geometries 

• Second-order central difference for both inviscid and viscous terms 

• Fully coupled, fully implicit algorithm 

• Efficient convergence for wide range of Mach numbers (from M < 10 10 to 
supersonic) 

• Finite-rate chemistry, realistic thermophysical properties 

• Multi-block and body-fitted curvilinear coordinates for complex geometries 

• k — c turbulence model 

• Stochastic liquid spray model (dilute spray), vortex model for droplet internal 
circulation and diffusion 

• Still a research code, require experience and knowledge of CFD and flow 
physics to use. 

Future Plans 

• Development of a more efficient solver 

• Extension to 3-D 1 

• PDF model for turbulence/chemical reaction closure 

• Thermal radiation model 

• Detailed soot and NOx kinetic models 

• Multi-grid and unstructured grid capabilities 

• Dense spray and high pressure (near-/super-critical) spray models 


Difficulties with Compressible Flow Algorithms at Low Mach 
Numbers 


• Disparities among system’s eigenvalues (stiffness), u, u + c, u — c, re- 
sulting in significant slowdown in convergence rate. 

• Singular behavior of pressure gradient term in momentum equations 
as Mach number approaches zero, 


4 *2 | 

p u + 


_P 1_ 
7M 2 


As Mach number is decreased, pressure variation (Ap* oc M 2 ) becomes 
of similar magnitude as roundoff error of the large pressure gradient 
term (p*/ 7 M 2 ). 



METHOD OF APPROACH 


Pressure Singularity Problem 

• Pressure decomposed into two parts: 

p=p 0 + Pg 

p g replaces p in momentum equations and retains p g as one of the 
unknowns. 

. Employs conservative form of governing equations, but uses primitive 

variables 

(p g ,u,v,h,Yi) 

as unknowns. Conservation property preserved and pressure field ac- 
curately resolved for all Mach numbers. 


Eigenvalue Stiffness Problem 

• Pressure rescaled so that all eigenvalues have the same order of mag- 
nitude Physical acoustic waves removed and replaced with pseudo- 
acoustic waves which travel at speed comparable to fluid convective 

velocity. 


EIGENVALUES RESCALING 


= U,U, i [l7(l + ± \Ju * (1 - ^) 2 + 4 ^(« i + °a) ] » u ’ u 


<*! = £*, a 2 — tyi U = a 1 u + a 2 v. 

For well-conditioned eigenvalues, scaling factor 0 taken to be 


0 = u 2 + v 2 . 


LIQUID-PHASE EQUATIONS 


Droplet Motion Equations 

(2) 

(3) 

(4) 

(5) 

Rep = 2 r -^[(u a -u p )' + (v,-v p n'\ 

Pa 

04 Rc 2 / 3 

(1 H ■¥—) for lie,, < 1000, 

iXCp t) 

C D = 

0.44 for Re p > 1000. 

Droplet Heat and Mass Transfer Equations 

^ = 2JV,/n(l + B), (6) 

pD f 

hdf _ 2NMI ± B) Lc ~ l . 

k ~ [(1 + B) Le ~' - 1] { ’ 


du p 

dt 

dv p 

dt 


dx p 

dt 

dy P 

dt 


= u 




= u 


p> 


3 C d fig Rc>p 


Pv r l 


16 

3 Cs£)fl>gRCp 


K - u p)» 


16 


Pv r l 


( V 3 - V p)» 


OMeRel^Pr 1 ' 3 

N, = 1 + c - 


[1 + 


1.232 ll/2 
Rc p Pr*/* i 


N p = l + 


0 .Zmtey'Sc'l* 

p+kSSp! 1 '” 


_ Yfgp — Yfg y X; Sp Wf 

" 1 -Yj„ i9P X }sp Wj + (l- X, SP )W a 



Droplet Internal Temperature Equations 
(Vortex Model) 


<Hp = 

dt r — 2 


’cfe l “^ +(l + c 


( 8 ) 


ar„ _ 

a = 0 ’ a^“i7 ( k )r * dt ’ 

i dT v _ 3 5I>. 

a “ *’ aa 16 rp dr )s 


NUMERICAL ALGORITHM 


Gas-Phase - ALLSPD code 
Liquid-Phase 

- Droplet motion equations (ODE) - Runge-Kutta method. 

— Droplet internal equations (PDE) - implicit method 
(Thomas algorithm). 

— Determination of spray time step for integration. 

- Stochastic separate flow model. 

Interaction Between Two Phases 


SPRAY TIME STEP 


• Droplet Velocity Relaxation time (< r ) 

tr = T<r)(?)( c ^r'- 

3 Pa v 

• Droplet Life Time ( t \ ) 


*1 = 


Zrh'Lpi 


• Droplet Surface Temperature Constraint Time (t,) 

I. = '"( 7-W/ A ' 

i r * 1 




• Local Grid Time Scale (i 5 ) 

• Turbulent Eddy-Droplet Interaction Time (t,-) 

= t e , z/ L e > t\u - Up | 
i; = if L e < t\u' - Up\ 

where 

U = Cj'^/e 
f, = L./(2<c/3) , «. 
i, - -r/n[l - i,/(r|u" - 0^|)] 

• Spray Time step - A T, pr 

Af 4 p r — O ^flm 




INTERACTION BETWEEN TWO PHASES 


1. Initialize gas and liquid phase variables. 

2. Solve liquid-phase equations. 

3. Evaluate spray source term, Hj. 

4. Solve gas-phase equations and update gas-phase 
variables, 

5. Update spray source term, Hj ? 

No, go to step 4. 

Yes, go to step 2. 


NUMERICAL TEST RESULTS 


Turbulent Backward- Facing Step Flow - non-reacting. 
Evaporating Turbulent Spray Flow. 

Gas Turbine Spray Combustion Flow. 



PARTIC1J: TRACES 


0.052 MACH 

g^90.00 DEG ALPHA 

tWxlOO GRID 











5.33 7.67 =.55 10.35 15.66 











OPERATING MODE 1 

NO ATOMIZING AIR 

LIQUID FLOW RATE : 1.26 grams /sec 
AIR FLOW RATE : 0.00 grams / sec 
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Convergence history for both non-combustion and spray combustion cases for the gas 
turbine spray combustion flow. 


1.2 norm 



CONCLUSION 


• ALLSPD code efficiently coupled with the SSF spray 
model. 

• Satisfactory convergence property for flows with and 
without spray. 
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Convergence history for EEE combustor 
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Chemically Reacting Flows 
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Part A. On Turbulent Transport of chemical Species in 
Compressible Reacting Flows 

Contents 

Mixing of chemical species in reacting flows 
Analysis 

Chemically reacting, turbulent flow equations 

Density-weighted, time-averaged Navier-Stokes equations 
Multiple-time-scale turbulence equations 
Species conservation equations for reacting flows 

Numerical Method 

Mixing and combustion of hydrogen in vitiated supersonic airstream. 

Comparison with measured data and other numerical results obtained using 
k-8 turbulence models and a PDF. 


Conclusions and Discussion 


Mixing of chemical species in reacting flows 


Chemically Reacting Laminar Flows 

Numerical (or semi-analytical) methods 


O.D.E. solvers 

1-D Numerical methods 

Numerical uncertainty is minimal 


Chemical kinetics 

Chemical kinetics (finite rate chemical kinetics & Reduced finite rate chemical 
kinetics) for certain fuels (hydrogen and some of hydro-carbons) have been 
tested and validated repeatedly. 

Numerical and theoretical analyses yield accurate results. 


Chemically Reacting Turbulent Flows 


Numerical methods 

Boundary layer or Navier-Stokes equations solvers 
(Except for semi-analytical analyses of PSR cases) 

Uncertainty is caused partly by numerical methods. 

Chemical kinetics 

Finite rate chemical kinetics & Reduced finite rate chemical are used. 

Uncertainty is caused partly by over-simplified chemical kinetics. 

Turbulence Equations 

Numerical analysis yields not so accurate results. Unceratinlty caused by turbulence 
equations is by far greater than that caused by chemical kinetics.* 

The probability density functions are used to improve the predictive capability of 
turbulent mixing of chemical species. However, some pdf methods yield physically 
incorrect numerical results. 

It is certainly important to better understand the turbulent mixing of chemical 
species in reacting flows. 

♦Westbrook and Dryer, Prog. Energy Combustion Science, vol. 10, pp. 1-57, 1984. 


Turbulent Mixing of Chemical Species in Reacting Flows 


1. Large eddy mixing 

2. Turbulent Mixing 

3. Molecular Mixing 

The use of accurate chemical kinetics alone can not yield accurate results for 
numerical analysis of turbulent reacting flows unless turbulent mixing is resolved 
correctly. 


Large Eddy Mixing 

Caused by separated flows, recirculating flows, and organized structures. 

The extent of calculated recirculation zone or the organized structures depends 
on the accuracy of both the numerical method and the turbulence equations since 
they are nonlinearly coupled with each other. 



Therefore, the capability to correctly resolve the large eddy mixing depends on 
the accuracy of the numerical method and the turbulence equations. 

Numerical investigations carried out during last decades show that k-e, ARSM, 
and RSM do not yield accurate results for complex turbulent flows. 

On the other hand, the multiple-time-scale turbulence equations* yield highly 
improved numerical results for various complex turbulent flows. 

*Kim and Benson, Int. J. Heat Mass Transfer, vol.35, pp. 2357-2365, 1992. 



Turbulent Mixing 


Caused mostly by energy-containing eddies and partly by fine-scale eddies. 

The capability to resolve the turbulent mixing of chemical species depends on the 
capability of the turbulence equations to correctly describe the turbulent 
transport of scalar variables (i.e., heat transfer, turbulent kinetic energy, 
concentrations, and convection-diffusion of species, etc) 

Need to be able to describe the chemical reaction-turbulence interaction (i.e., 
turbulent mixing is enhanced and shear layer thickness is widened by chemical 
reaction)* 

The multiple-time-scale turbulence equations can resolve the cascade of 
turbulent kinetic and the nonequilibrium turbulence phenomena. 

(Present numerical results show that the M-S turbulence equations can resolve 
the chemical reaction-turbulence interaction.) 


* C.T. Chang et al., "Comparison of reacting and non-reacting shear layers at a high 
subsonic mach number," AIAA paper 93-2381, 1993. 


Molecular Mixing 
Caused by molecular diffusivity. 

Theoretical and numerical analyses of chemically reacting laminar flows yield 
accurate results. 

Various molecular diffusion equations can accurately describe the molecular 
diffusion of species. 

The Lennard- Jones 12-6 potential law is used in the present study. 

Various numerical methods that yield accurate results for laminar flows can 
resolve the molecular mixing. 


Multiple-time-scale turbulence equations 
Energy containing eddies 
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Turbulent Eddy Viscosity 

p t =pcpk 2 /£p where c^ = c pf £t/£ p 

Remark: Single-time-scale turbulence models can not resolve nonequilibrium 
turbulence phenomena. 


Species conservation equations for reacting flow's 
Chemical species concentration equation 

i(pYi) + V- {pY i(v + V,)}= Wj 

where the diffusion velocity, Vi, is approximated using the Fick's law given as 
Y i V i = -(D iJ + D i , t )VY i 

The production rate of the i-th species, Wj, is given as 

. N r , , 

wi= I MiJn^-n^Jwk 
k=l 

where 

Ns v/ N S 

(Ok=^f,k n c# - ^b,k n c > k 

j=l 3 j-1 3 



Chemical reactions for the combustion of H/j in a vitiated supersonic airstream 
are described using 9 chemical species (H 2 , 0 2 , H 2 0, OH, O, H, H0 2 , H 2 0 2 , and 
N 2 ) and 24 pairs of reaction-steps (Burks and Oran, 1981; Kumar, 1989). 

A fast chemistry can not be used to describe the fine details of chemically 
reacting flows. 

A reduced chemical kinetics can not be used confidently due to the uncertainty 
contained in the reaction mechanisms. 

The use of a detailed finite rate chemistry may make it difficult to obtain a fully 
converged solution due to the coupling between the large number of flow, 
turbulence, and chemical equations. The numerical method needs to be 
strongly convergent. Accuracy also depends on the capability of turbulence 
equations used. 


Numerical Method 

The numerical method is a finite volume method that incorporates a 
pressure-staggered mesh and an incremental pressure equation for the 
conservation of mass. 

Predictor Step: Solve momentum equation. 

(pCl+A*)u**=X a£u£% S*- + PC2U]! 1 ' 1 - PC3UP* 2 (1) 
nb ax i 

i 

Corrector Step: Correct the velocity field to be divergence free 
Incremental pressure equation 

dip' d I 1 9p’|_ ^ u j 

c)xj\RT -i ) 3xj|(pCi+A*)9xj| dxj 

Incremental velocity equation 

U ' = 1 ^ 

1 (pQ+A*) dxj 


(3) 



Velocity and pressure corrections 


ur=ur+u'i 

(4) 

p** = p* + p' 

(5) 


Solve eqs. (1-5) iteratively until all flow variables are converged. 



Combustion of H 2 in vitiated supersonic airstream (Burrows and Kurkov, 1973) 





- RPLUS (w/o PDF) 





Total temperature 



Calculated OH contour. 

The numerical result obtained using the M-S turbulence equations indicates 
that the combustion occurs at approximately middle of the channel. The 
calculated flame location is in excellent agreement with that observed in the 
experiment. 






Temperature contours obtained using the M-S turbulence equations 


It can be seen in the temperature contours that the calculated flame location 
obtained using the M-S turbulence equations are in correct agreement with the 
measured data. The contours also indicate that the temperature increases 
within in a short distance. The trend is in correct agreement with experimental 
observations that temperature increase occurs within a finite flame thickness. 



Temperature contours obtained using the pdf* 


The pdf fails to predict the correct flame location (i.e., ignition delay). The 
slowly increasing temperature field indicates that the pdf may not be able to 
predict a correct flame front. The numerical results obtained using the pdf are 
not in correct agreement with the physics of combustion. 

* A. T. Hsu, Y.-L.P. Tsai, and M. S. Raju, "A PDF approach for compressible 
turbulent reacting flows/ 1 AIAA Paper 93-0087, 1993. 



Conclusions and Discussion 

The calculated species concentration profiles are in as good agreement with 
the measured data as those obtained using the pdf. 

The flame location (ignition delay) obtained using the multiple-time-scale 
turb ulence equations is in excellent agreement with the experimentally 
observed onset of ultraviolet radiation. 

Cascade of the turbulence filed is influenced by the extra strains caused by 
chemical reaction. 

Both the numerical results and the measured data exhibit enhanced mixing 
of the hydrogen and vitiated airstream for the reacting case. 

The M-S turbulence equations can resolve the chemical reaction-turbulence 
interaction. 

The pdf produces slowly increasing temperature field and it fails to predict 
the ignition delay. Thus the numerical results obtained using the pdf are not in 
correct agreement with physics of combustion. 


Part B. Unsteady Transitional Flows over Forced Oscillatory Surfaces* 

Contents 

Nomenclature 

Unsteady turbulent flow equations for flows with moving boundaries 

Navier-Stokes equations defined on Lagrangian-Eulerian coordinates 
Multiple-time-scale turbulence equations 

Numerical results 

Unsteady transitional flow field and comparison with measured data. 


* S.-W. Kim, K.B.M.Q. Zaman, and J. Panda, "Calculation of unsteady transitional 
flow over oscillating airfoil" in Separated Flows, eds. J.C. Dutton and L.P. Purtell, 
Proceeding of ASME Fluid Engineering Conference, Washington D.C., June 20-24, 
1993. 



Nomenclature 



Lagrangian-Eulerian coordinates 

x: fixed reference coordinates 
moving coordinates 
J = id^j/dajl 


Unsteady transitional flow equations with moving boundaries 


Conservation of mass equation 

d(_r\ T 3 


Conservation of linear momentum equation 

|( pu i J ) = J jr(p u i( u f - U j)} + J a^‘ J || 

Convection-diffusion equation for scalar variables (i.e., <|) = |kp, £p, kp £[, etc. | 


o 




a<t> 


■W) 


Numerical Method 


The unsteady transitional flow equations are solved using the same finite 
volume method. Time-integration is made using an iterative-time-advancing 
scheme. 

Comparison of Unsteady Flow Solution Techniques* 

Iterative Time-Advancing Scheme (ITA). 

Simplified Marker and Cell (SMAC). 

Pressure-Implicit Splitting of Operators (PISO). 

The ITA that can best resolve the nonlinearity of the Navier-Stokes equations, 
and yields the most accurate results. 

The SMAC is the most efficient computationally and yields accurate numerical 
results for laminar flows. 

The PISO is the most unstable numericaly and yields less accurate results. 


* Kim & Benson, Computers and Fluids, vol. 21, pp. 435-454, 1992 



(a) a =15°/ t/T = 0.0 


(b) a = 25°' t/T = n/2 


(c)a=5°/ t/T = 3/2 n 


Oscillating airfoil and moving mesh 




t/T=0.0 t/T=0.2 
(5°u) (ll°u) 


1234567 5678 5678 5678 5678 5678 

u (m/sec) 

Ensenble-averaged velocity profiles at x/ c=1.0 
The deteriorated comparison at a =19°d 

(i) The hot wire can not accurately measure the velocity components when the 
flow is misaligned more than approximately 30° from the hot wire axis. 

(ii) The interaction between the DSV and the TEV occurs in a relatively coarse 
mesh region and the numerical method yields somewhat deteriorated results 











Turbulent viscosity contour ( = 10) 



CONCLUSIONS AND DISCUSSION 

Numerical method successfully predicts the Dynamic Stall Vortex and the 
Trailing Edge Vortex. 

The calculated ensemble-averaged velocity profiles are in good agreement 
with the measured data. 

Both the numerical results and the measured data show that the transition 
from laminar to turbulent state and relaminarization occur widely in space and 
in time. 

The good comparison between the numerical results and the experimental 
data are attributed to the capability of 

(i) the ITA that can best resolve the nonlinearity of the Navier-Stokes 
equations, 

(ii) the new pressure correction algorithm that can strongly enforce the 
conservation of mass, and 

(iii) the Multiple-time-scale turbulence equations that can resolve the 
transitional nonequilibrium turbulence field. 
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FOURIER ANALYSIS 


Fig. 2 Dispersion characteristics of the CD6-RK46 scheme for CFL=0.6. 







A LAMB VORTEX IN AN UNIFORM FLOW 



wave mode 


SHOCKED SOUND WAVE 



with A I 






Fig. J | Tlic vorticity magnitude contours for the vortex pairing. 










The Fourier coefficients of the instability waves in a compressible 
layer perturbed by the most unstable mode and its subharmonic. 




CONCLUDING REMARKS 
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Turbulent Back-Facing Step 
Flow and the k-e model: 

A Critical Comparison 




ref: A\ AA-iynz? 


Christopher J. Steffen, Jr. 

Center for Modeling of Turbulence and Transition (CMOTT) 
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INTRODUCTION 

CMOTT at NASA Lewis: 

to address propulsion related turbulence issues 

□ assessment of state-of-the-art models 

□ development of new models 

□ validation of new techniques 



Objective of this work: 

to assess the performance of several low Reynolds number k-e 
formulations compared to the standard high Reynolds number 
form for separated flow over a step 


Lewis Research Center 


CFD Branch 



Need to clarify this issue for code developers 
Critical comparisons between k-e models for: 

□ Channel flow: Lang & Shih (1991) 

□ Flow past a hill: Michelassi & Shih (1991) 


Awa (1990) compared the high Re t and Chien k-e models for 
several flows: noted the deficiency of Chien for step flow skin 
friction results . . . 


Contradicts idea of “more work<~*-better result” 

Is this a property of all Low Reynolds Number k-e models? 
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k-e TURBULENCE MODEL 

Specifically, v t = 

Model transport of k and e: 


fc, t+Uj&j- 

e ,t+Uj€ J '- 


v+— 1 * i 




= V—e 4- D 


j j 


, ^ i 


J J 


= CJi—V — 


where the production of k is: V= — (uiUj)|(Uij 4- Uj,i). 


Primary difference between formulations: 
definition of / m ,/i,/ 2 ,D,E 


Lewis Research Center CFD Branch 

Low Re t forms: 

□ numerical domain extends to wall 

□ near wall profiles of U, k, and e are resolved not assumed 

□ typically first cell located y + « 1.0 

□ damping functions (fj, f 2 , and f^) and additional terms (D 
and E) tuned to specific flows 

Generality suffers somewhat . . . 

Formulation of above terms can affect the generality as well: 

□ JL and LS are functions of dependent variables alone 

□ CH and SL are functions of both dependent and independent 
variables 
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BACK-FACING STEP (BFS) FLOW 

Experiment of Driver and Seegmiller (AIAA J., ’85) 
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Configuration: 

□ Re st ep of 33420 

□ inlet Mach Number of 0.128 

□ inlet tunnel 80H 

□ exit tunnel 60H 

Chosen for: 

□ LDV data for velocity, turbulent stresses 

□ wall static C p and wall Cf 

□ reattachment length (time averaged) 


yf 



NUMERICAL EXPERIMENT 

Code: modified version of DTNS2D (Gorski ’88) 


Method: 

□ pseudo-compressibility (Chorin ’67) 

□ upwind differenced, TVD scheme for convective terms 
(Chakravarthy & Osher ’85) 

□ approximate factorization for time integration 

□ multiblock configuration 

□ decoupled (lagged) treatment of turbulence equations 

□ implicit treatment of turbulent source terms 


Lewis Research Center 


CFD Branch 
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Velocity profiles: 


o Driver and Seegmiller 

Jones and Launder 

Launder and Sharma 

Chien 

Shih and Lumley 

Standard High Re, 
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Wall static pressure coefficient: 


o Driver and Seegmiller, stepside wall 
a Driver and Seegmiller, opposite wall 

Jones and Launder 

— Launder and Sharma 

Chien 

Shih and Lumley 

Standard High Re, 
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Turbulent kinetic energy profiles; 


o Driver and Seegmiller 

Jones and Launder 

Launder and Sharma 

Chien 

Shih and Lumley 

Standard High Re, 
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Turbulent shear stress profiles; 


o Driver and Seegmiller 

Jones and Launder 

Launder and Sharma 

Chien 

Shih and Lumley 

Standard High Re, __ 
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Primary reattachment length. x r : 


[model 

JL 

LS 

CH 

SL 

HR 

\x r 

4.9 

5.4 

5.4 

5.1 

5.5 

■— ; ~ i = 


Driver & Seegmiller recorded: 6.26H 


Secondary reattachment length. x r : 


[model 

JL 

LS 

CH 

SL 

\ x r2 

0.0 

0.0 

0.8 

0.1 


Driver & Seegmiller recorded: 0.5 < x r 2 <1.8 
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CLOSING COMMENTS 


□ JL, LS, CH, SL and HR k-e models tested for back-step flow 

□ In general, results very similar between 5 different forms 

□ All do reasonably well for U and C p , albeit LR slightly 
better 

□ Turbulence profiles possess misplaced peak values as well 
as overpredictions in recirculation zone, common to all 

□ Reattachment length and C t are difficult to properly resolve 
(14% error) especially for LR forms 

□ Definition of secondary recirculation zone inconclusive on 

this grid 

□ Perhaps a sensitivity to pressure gradient would help the 
LR damping functions . . • 
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SUBSONIC INLET FLOWS WITH TRANSITION 


BY 


Danny Hwang and Kyung Ahn 


September 16, 1993 
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INTRODUCTION 

Real flow simulation is needed 

(laminar region + turbulent region) 


for: 

(1) cruise drag prediction 

(2) separation angle of attack prediction 



Difficulties: 


• Transition is a function of many variables 

- free stream turbulence 
-surface roughness 
-streamwise surface curvature 

- pressure gradient 
-etc, 

• Skin friction coefficient which is the best 

transition indicator is difficult to calculate 


METHODS OF TRANSITION PREDICTION 

• Euler/boundary layer approach 

• Euler/boundary layer/stability analysis 

• Navier-Stokes calculation with turbulence models 

- Baldwin-Lomax with transition 

- RNG turbulence model 






TRANSITION PREDICTION 


• Flat plate 

• NACA0012 


• ADP inlet 


pnickpiotec.q, BL w/tron model 













LAM 1 N A R/T l J R RULE N 'I' T RA N SI TI ON 




Laminar/ i u rbulent Transition 


Laminar region I urbulcnt region 



Laminar region prediction: 

Lewis’ s boundary layer analysis : 8.17 % of chord length 
Langley’s stability analysis : 8.00 % of chord length 

On going projects: 

RN(« turbulence model in I’AIO 3D 
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Several methods have been tested to predict the 
transition of an ADP inlet. The prediction of the 
transition at cruise conditions agrees very well 
among the methods investigated except for the 
Baldwin-Lomax/Transition model which predicts 
early transition. The accuracy of prediction 
methods will be judged by comparing with 
upcoming laminar flow experiment. 
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Analysis of Supersonic Flows using k-£ Model and the 
RPLUS code; Progress towards High Speed Combustor 

Analysis. 


by J. Lee 
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for Workshop on Computational Turbulence Modeling 
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Problem of Interest - High Speed Combustor Flow Fields 
Parameters need to be Resolved 
Key Problems of Interest 
k-e Model and RPLUS code 
Numerical Technique 
Models being Tested 
Some Results 
Summary 


Sverdrup 

Technology Inc. 
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Problem of Interest 



Analysis of Chemically Reacting flow inside of Supersonic RAM jet 
Combustors-Two Key Parameters need to be determined. 
Mixing/Combustion Efficiency 
Kinetic Energy Efficiency (Flow Losses) 

Inlet, Diffuser, etc.. 

In order to do get some ideas on those parameter following (Potential 
Loss Mechanisms) must be modeled/determined correctly. 

Mixing, Shear, 

Turbulence, Vorticity, 

Shock-waves, Heat Transfer, 

Fuel Injector Drag, Poor Wall Pressure Integral, 

Chemical Dissociation. 

from 2nd JANNAF workshop on SCRAMjet Combustor performance workshop 
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Mixing and Injector Design 

At High Mach Number(M ~ 5.0 +). 

Doesn’t mix well! 

The Natural diffusion mechanism very INEFFECTIVE. 

Fuel Residence Time Extremely Small- Even with Fast Fuel Such as H 2 

. Geometrical Complexities x 

To induce Favorable mixing and Flame holding features 

Back-Step/Stream Wise Vorticity/Shock-Wave Interactions 
Unsteady Mechanism also being Envisioned as mixing enhancement 
Kumar, Bushnell and Hussani(1987) 




Sverdrup . . 

Technology Inc. Internal Fluid Mechanics Division 


Introduction of Externally Generated Mixing Enhancements 

Some External helping hand needed => Modeling Difficulties. 
Externally Generated Vorticity Through Sweep angle of the Ramp 
injector. 

Davis(l990), Riggins and McClinton(1990), Drummond(1991). 

. Multiple Transverse Injection. 

Hartfield et. al. (1991) 

. Flame holding tricks/ Back-step with Recirculation. 

Hartfield et. al.(1991) . 

. Simplified analysis of these features very difficult because of limited 
database/understanding (Attempts are being made using CtD 
solutions- JANNAF Combustor Subcommittee). 
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Numerical Modeling(CFD) of Combustor Flow Field 


• CFD Analysis. 

• Numerical Modeling=> Overall Analysis of performance => Difficult 

• Overall Laminar Flow Fields with Complex Geometry/Finite Rate 
Chemistry has been demonstrated. 

• Finite Rate Chemistry Model^ Yoon and Shuen(1989) 

• Multiple Grid Blocks- Moon (1991) 

• Analysis of a typical Injector Configuration with Zero Equation 
Turbulence Model using LU Scheme(RPLUS) code- Lee(1993) 
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Simple Zero Equation Turbulence model with multiple wall scaling 
Buleev-Inverse square rule can be used to extend model in to three- 
dimensional form. (Lee (1993)) 

• Good News/Bad News 

• Typical velocity profiles can be reasonably predicted. 

• Over all combustor flow features can be reasonably predicted. 

• Near-wall temperature characteristics near non-equilibrium region 
around the injector and separated flow were poorly predicted. 

• Overall spreading behavior of shear region poorly predicted. 

• Two Equation Transport Turbulence Model has the potential to ease 
some of these difficulties. 
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Lewis Research Center 


1st bow Shod/ / 2nd Bow Shock 


Normolized V-velocity 


Normoli7ed IJ- velocity 

— RPLUS(IW) 

_ - - RPLUS(2W) 

§5 RPLJSC3W) 

i ODotc McDonicI 

1 .4 
1? 

... RPLUS(IW) 

RPLUS(2W) 

-- RPIUS(3W) 
QOoto McDaniel 



Normolired Temperoture 


- rpluS(iw) 
- RPLUS(2W) 

- RP| US(3W) 



Normolized Pressure 

— RPiusnw) 

RPLUS(JW) 

RPtUS(3W) 

QUoto- WcDonie! 



Sygrd rup Internal Fluid Mechanics Division 


Lewis Research Center 
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ye.. LAi. 



Two Equations Transport Tlirbulence model are being Analyzed 


• High speed turbulence models are some what Deficient (The 
deficiencies are well documented(Marvin(1986), Wilcox(1993)). 

Effect of Compressibility 
An-isotropy (Low/High Speed). 

Non-Equilibrium Flow Features (Low/High Speed). 

Near- Wall Flow(Low-Reynolds Number Features (Low/High Speed)). 
Inflexibility of handling Complex Geometry- Invariance Principle 
(Low/High Speed) 

Large Dependence in the Numerical Methods Used 
(especially elliptic Solvers). 

Appropriate Initial/Boundary Conditions 
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K-e Model-RPLUS Development 



LU Based k-£ Model Solver-De-coupled Approach. 
Mean-Turbulence Transport Equations 

LU-SSOR- Yoon and Shuen- Explicit Terms Centrally Differenced 
LU-SW -Steger and Wanning- Explicit Terms Upwind Differenced 
k-e Models 

Convective Terms + Diffusive Terms + Source Terms = 0.0 
Model Only differ in Low-Reynolds Number Character. 

Models performance are being Evaluated. 

Implicit Source Term Handling Strategy also Being Studied 
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Internal Fluid Mechanics Division 


k-£ Tbrbulence Models being studied for potential used in Three 
Dimensional RPLUS Code. 

• Low-Reynolds Number Model plus Dilatational Terms 

Chien (1976) 

Launder-Shima( 1 976) 

Shih(1990) 

Various CMOTT derivatives of k-e Model 
Realizability 
Invariance 

Simplified Boundary-Conditions 

• Performance of the Low-Reynolds number K-e model in low-Mach 
number flows have been demonstrated (Patel, Rodi and 
Scheuerer(1985), Steffen(1993), Launder(1992)). 

• Some of the Potential Difficulties in high speed turbulence model are 
well documented (Marvin(1993), Coakley and Huang(1992). 


it 
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Evaluation and Development of the RPLUS/k-8 Model Solver 

Various 2D-3D problems are being studied to optimize the numerical 
method and to Evaluate model performance in supersonic flows in 
context to the LU based numerical Technique. 

Simple 2D k-e models are also being used to study various 
components of the flowfield generated by the complex combustor 
geometry previously shown. 

Studying the Numerical method/Model Behavior/Model Performance. 
2D Supersonic Turbulent Boundary-Layer- Skin Fraction/Heat transfer 
(NASA Ames Database). 

2D Supersonic Shock- Wave Boundary-Layer Interaction- Skin fraction/Heat 
-Transfer/ Shock- wave( A. Smits (1990’s)) 

2D Shear-Layer - Mixing (H. Lai(1993)) 

3D Fin/Flat Plate Interaction- 3D Comer Flows-Interaction Developed through 
a Fin generated Shock- Waves. (D. Davis(1992)) 
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Supersonic Mixing Layer 

RPLUS vs. Dutton (case 2) 


. ^ 

- RPLUS 2 

C Dutton, case 2 * 

O 



(U-U2)/dU 
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Boundary Layer 
Mach 2.87 
Re/m = 6.3x1 0 7 /m 


Law of Wall Profile 


50.0 


40.0 


0 30.0 


LU-SSOR/Chieu 
LU *SSOR/L*urMk r-Sh»rma 
LU-SSO R/ShQi-Lumlcr 
Lee«C ftL 

DaU SmiU rt. *L 
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Validation of a k-e Model in RPLUS2D Code 
for Non-reacting/reacting Subsonic Shear layers 


H. T. Lai 

Workshop on Computational Turbulence Modeling 
September 15-16,1993 


OVERVIEW 

- To simulate an expirement at NASA Lewis by 
Marek and co-workers for subsonic shear layers 

- Computations for nonreacting/reacting flows 

- Grid refinement study for nonreacting case 

- Compressibility effects in k-e models (nonreacting) 

- Torch effects to sustain combustion (reacting) 



Sverdrup 

Technology Inc. 


Internal Fluid Mechanics Division 


Turbulent Shock- Wave/Boundary Interactions 
Mach 2.87 

Ramp Angle = 8.0 degrees 
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Other Factors 



Optimum Numerical Strategy with in LU frame work. 

Effects of Initial condition. 

Modeling of Compressibility terms/Dilatational terms. 

Modeling of Turbulent terms in the Finite Rate Chemistry Model. 
Anisotropy of Turbulence 

Effects Upstream and Down stream Influences (Inlet(K. Kapoor) and 
Diffuser(?)). 

Chemistry-Turbulence Model Interactions (A. Hsu-PDF). 

Numerical Robustness(A. Suresh). 





The RPLUS Code 


- Navier-Stokes and species conservation equations 

- Finite— volume scheme 

- Options for central or upwind differences 

- Jameson— type dissipation 

- LU- ADI algorithm 

- Vectorization on oblique grid lines/planes 

- Hydrogen-air chemistry 

- implicit chemical source terms 

- High Reynolds number k-e model 

- k-e equations uncoupled from other equations 


NONREACTING SHEAR LAYER 


U=140 m/s, T=304.3°K, P=1.013xl0 5 N/m 2 

y=5.08cm 


ti=2.25%, AIR 

symmetry 

U=390 m/s, T=822° K, P=1.1013xl0 5 N/m 2 
ti=2.25%, AIR 

y=-5.08cm 

„ cr 


x=55.9cm 


GEOMETRY AND FLOW CONDITIONS 



GRID DISTRIBUTION 


NONREACTING SHEAR LAYERS 


o MEASURED 

RPLUS 

WEC 



U. m/w 


PROFILES OF AXIAL MEAN VELOCITY 



PROFILES OF TOTAL RMS FLUCTUATION 


NONREACTING SHEAR LAYER 



M=.679 


MACH NUMBER CONTOURS 



TURBU LEN T KINETIC ENERGY CONTOURS 



DISSIPATION RATE CONTOURS 






o TV K 2 3 

DISSIPATION RATE 3X10 ° m 

PROFILES OF DISSIPATION RATE 



SHEAR LAYER THICKESS VARIATION 



GROWTH RATE VS. COMPRESSIBILITY EFFECT 



REACTING SHEAR LAYER 


U= 135 m/s, T=366.5°K, P=1.013xl0 s N/m 2 

— 

ti=5%, 3.94% H 2 and 96.06% N 2 


U=385m/s, T=810.9°K, P=l.Ol3xlO s N/m 3 
ti=6%, AIR 



FLOW CONDITIONS 


^ dimensions 1cm 

hydrogen torch ti=5%,T=1250° K, 1.3% H 2 + 21% 0 2 + 77.7% N 2 
V=0 (parallel inflow) 

pARALLEL jnplqw 


hydrogen torch 


dimension=2. 12cm 

ti=l 1%, T=1250°K, 0.1% H 2 - 

V=50m/s (oblique inflow) 


• 21% 0 2 + 78.9% N 2 


OBLIQUE INFLOW 


REACTING SHEAR LAYERS 

O MEASURED 

RPLUS(OBUQUE) 

RPLUS(PAJFULLEL) 

WE C- PDF 



PROFILES OF TOTAL RMS FLUCTUATION 









REACTING SHEAR LAYERS 
TEMPERATURE CONTOURS 



T=300° K 

T=1SKX) K 


PARALLEL INFLOW 



OBLIQUE INFLOW 


REACTING SHEAR LAYERS 
H 2 0 MASS FRACTION CONTOURS 



OBLIQUE INFLOW 




CONCLUSIONS - 

- Growth rate underpredicted, especially reacting 

- Large decay of k near shear layer edges 

- Small production of k or large dissipation rate e 

- But center peaks in k are constant and overpredicted 

- Small improvement with grid refinement 

- Inhibition of growth with compressibility effect 

- Significant torch effects for reacting flows 


0 100 200 300 „ 400 500 600 

X, mm 

VARIATION OF WATER VAPOR MASS FLOW RATE 


Computational Aero*Acoustics as a Branch of Turbulence Research 
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• The initial region of jet is dominated by large-scale, wave-like coherent structure, which 
is believed to be the dominant sound source. 





Ficcrr 13. Maximum large-scale Unicture amplification aa a function of Strouhal number; companion with 
Moore** ( 1977 ) experimental data. 


Calculated Spectra of Sound Intensity in Decibels Referred to 
10-12 w m -2 Due i 0 Coherent Structures 
at Various Emission Angles for U e = 1 95 m s* 1 
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The Study Indicates: 

II. Large-Eddy Simulations For Prediction of The Sound Source in a Supersonic Jet ^ Jhc |arge . sca | e structure seems to be the dominant sound source. 

• Results are sensitive to approximation in the sound source. 

<• cmirrr— But. the resolution requirements • Lighthill’s theory predicts some results consistent with observations and some are not. 
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-Real pari of u at several Slroohal number* under harmonic Inpot disturbances. 





HI. Far-Ficld Sound 

(1) Extend Navier-Stokes computational domain to (lie far Held: 

(a) Prohibitive storage requirements 

(b) acoustic scales are different from fluid scales 

Strouha 

(2) LighthilPs acoustic analogy mimhnr 



(3) Finite-difference approach — Linearized Euler equations 
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(4) Kirchhoffs method 



• Predicted directivity of jet noise seems to be consistent with observation. 


JET NOISE PREDICTION USING A k-e TURBULENCE 

MODEL 0W 


A. Khavaran 

Sverdrup Technology, Inc. 
NASA Lewis Research Center 


Lewis Internal Workshop 
on 

Computational Turbulence Modeling 
September 15-16, 1993 


MODELING APPROACH 


Source Spectrum Calculations 

- Acoustic Analogy 

- Ribner and Batchelor Assumptions 

- Calculation of the Source Spectrum and its Characteristic Frequency Based 
on Time-Averaged Flow Calculation ( with a k-e Turbulence Model) 


Sound/Flow Interaction 

- High Frequency Asymptotic Solution to Lilley’s Eq. for Multipole Sources 
Convecting in an Axisymmetric Parallel Flow ( Balsa & Mani) 



GOVERNING EQUATIONS 


Lighthill’s Equation 


— - c 2 V 2 p = 

dt 2 P dXidXj 


Tij = pViVj + 6jj(p - c 2 p) - e {j 

,dV.dV i _2 cm 

eij ~^dxj + dxi 3 6,j dxj 

— 0(4"). fie= — 

pViVj 'Re P 

Idp-dp = (^)* 


The effects of source convection and refraction are included in the source term 


Lilley’s Equation 


D D 2 a 
DV Dt 2 


d 2 dfx 2 — 

dxi° dx, dxidxj ° dxi 


2 dV i dV !i dV L + 

dxi dxj dxk 


D . D 1 Ds 

— I— — ( — — )] + viscous terms 

Dt l Dt'c v Dt n 
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7 Po 


D_ _ d_ d 
Dt dt k dxk 


• The effects of source convection and refraction are included in the operator term of 
Lilley’s eq. 



g^.w-rg Spec trum Calculations 


d pressure autocorrelation in the far field due to a fimt 
Mean-square sound pressure shielding) 

turbulence (in absence of convert, on and flu,d sh,eld g) 



Fourth-order velocity correlation tensor 


/ +oo 

{ViVj){v' k v'i)dt 
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Source strength (Quasi-incompressible turbulence) 



Reduction in order of correlation tensor 


Sijkt = SikSjt + SitSjk + Si j Ski 



Separable second-order tensors 


Sn(r,() = R tj (()G(r) 


Isotropic turbulence model of Batchelor 

= Te-*W L *f x {[1 - x(t/L z ) 2 }6 tj + *&/ L \ } 

t = | m, ( 2 = ti+ei+(i 


Gaussian correlation time delay 


G{t) = e -< r / T ">’ 


Source spectrum component 


W«)~ P 2 ^(ClT 0 ) 4 e = ia ^ 1 


L x ~ 
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2 V ' Vi 


Characteristic time delay of correlation 
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or r 0 
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Doppler shifted frequency 

n = 2 jt/C, C = ^(1 - ~McCosfff + 
M e = .5M + /3 c Mj 


Sound /F] ow Interaction 


Mean square pressure in the far field 

pi(R,d,Q.) = f A(a xx + 4a xy + 2a yy + 2a yz )dy 
Jy 


Source term 

(^fm 

A ~ - McosOf( 1 - MgCosdY 


Shielding function 


$*( r ) = 


(1 - Mco ^) 2 (%) 2 - cosr 6 

(1 - Mc.cosO)' 


M(r) = U{r)/Coo M c = UjCoo 


Design Pannatcrs for C-D No///lc 



Throat Diameter . 5.1 in. 

Exit Diameter 5.395 in. 

Distance' from Throat to Exit 5.525 in. 
Exit Velocity 2409 fps 

Ami > if li l Velocity 400 tps 

Pressure Ratio 3.121 
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7 . Turto. Intensity 


Comparison of the Turbulent Intensity Profiles With Data 
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OASPL DIRECTIVITY (ARC=40 ft) 



Angle From Jet Axis 



Design Parameters for C*D Nozzle 



Exit Diameter 5.395 in. 

Throat Diameter 5.1 in. 

Exit Vejocity 2409 fps 

Pressure Ratio 3.121 

Stagnation Temperature 1716 R 


NOISE SPECTRA ( CD Nozzle) 



Freq. Band No. Freq.BandNo. 


Band 15 is 125 Hz 






GRID 
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OASPL DIRECTIVITY (ARC=10 ft) 


125 



Pressure Ratio 3.61 

Stagnation Temperature 532 R 
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NOISE SPECTRA (CPN) 



Freq. Band No. 

Band 15 is 125 Hz 


6 = 60 ° 




SUMMARY 


Source strength was evaluated using the PARC code with a k-e turbulence model 

The characteristic Strouhal no. was obtained from k and e ( with an empirical 
proportionality constant) 

Time-averaged velocity and temperature predictions were used to assess the 
sound/flow interaction 

Constants used in the supersonic convection factor are determined empirically 

The SPL directivity and spectra demonstrate favorable agreement with data ( 
specially at angles close to the jet axis) 

The empirical constants used in these predictions need to be investigated for other 
geometries 






OUTLINE 


NUMERICAL SIMULATION OF UNSTEADY SUPERSONIC FLOW USING 
AN IMPLICIT ALGORITHM FOR THE STRONGLY COUPLED 
NAVIER-STOKES AND K-E EQUATIONS 

by 


S.H. Shih* A. Hamed # and J.J. Yeuan# 

* ICOMP , NASA Lewis Research Center 

# Dept, of Aero. Eng. & Eng. Mechanics 
University of Cincinnati, Cincinnati, Ohio 
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Fig. 41 Computational Grid for Supersonic Flow Over an Open Cavity. 
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Schematic of the Solution Domain for Supersonic Flow Over an Open Cavity. 
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IMPLEMENTATION OF A k - 8 MODEL 
IN SPECTRAL ELEMENT CODE 


Wai-Ming To 


Sverdrup Technology, Inc. 


OBJECTIVES 


• Study the flow physcis in problems that are generic to 
turbomachinery flows. 

• Assess accuracy of turbulence models. 


EQUATIONS OF MOTION 



TIME DISCRETIZATION 

• 2nd order predictor-corrector 

• Capacitance matrix method to solve for pressure 


TIME DISCRETIZATION (cont.) 


u 


( 1 ) _ 


Ui +Atfi-At^^-+At 3 


dx; 


dx. 


du) 


( 1 ) 


dx; 


= 0 




du? a«7 
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dx; dx; 
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TIME DISCRETIZATION (cont.) 
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where p = + p n 



TURBULENCE MODEL 


Eddy viscosity concept 

subgrid model — » 2-equation model 

Zero equation model 


v, = K 2 y 2 1 - e y 726 


v t = 0.018 Moo 5* 


TURBULENCE MODEL (cont.) 


Wall boundary conditions 


K>? —In Ey + 
K 


T + hvp) 

where y = - 


K V/ 
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A Comparative Study of Turbulence Models in Predicting 
Hypersonic Inlet Flows 

by 


Kamlesh Kapoor 


TEST CASE AND TURBULENCE MODELS CONSIDERED 


• THE NASA P8 INLET, WHICH REPRESENTS CRUISE CONDITION OF 
OF A TYPICAL HYPERSONIC AIR-BREATHING VEHICLE. WAS 
SELECTED AS A TEST CASE FOR PRESENT STUDY. 

• PARC2D CODE. WHICH SOLVES THE FULL TWO-DIMENSIONAL 
REYNOLDS-AVERAGED NAVIER- STOKES EQUATIONS. WAS USED 
FOR THIS STUDY. 

• THE RESULTS ARE PRESENTED FOR A TOTAL OF SIX VERSIONS OF 
ZERO- AND TWO-EQUATION TURBULENCE MODELS. 

• ZERO-EQUATION MODELS 

• THE BALDWIN-LOMAX MODEL 

• THE THOMAS MODEL 

• A COMBINATION OF THE B.L./THOMAS MODEL 

• TWO-EQUATION MODELS 

• LOW-REYNOLDS NUMBER MODELS 

• THECHIEN MODEL 
•THE SPEZIALE MODEL 

• HIGH-REYNOLDS NUMBER MODEL 

* THE LAUNDER AND SPALDING MODEL 


EXPERIMENTAL BACKGROUND 

• THE EXPERIMENTAL INVESTIGATION OF THE P8 INLET WAS CONDUCTED 
AT NASA AMES' 3.5-FOOT HYPERSONIC WIND TUNNEL. 

• THE INLET WAS A MACH 7.4 RECTANGULAR MIXED COMPRESSION (WITH INTERNAL 
COMPRESSION RATIO OF 8) DESIGN WITH EXITING SUPERSONIC FLOW. 

• INLET COWL HEIGHT - 18.33 CM. 

•OVERALL LENGTH - 136.2 CM. 

• TEST CONDITIONS: 

•MACH NO -7.4 

•TOTAL PRESSURE - 4.14 x 10 6 N/nr? 

•TOTAL TEMPERATURE -81 1°K 
•REYNOLDS NO - 8.86 x 10 6 /m 

• MODEL WAS WATER COOLED AND ISOTHERMAL WALL CONDITIONS 
WERE MAINTAINED; THE WALLS TEMPERATURE - 302°K 

• THE TRANSITION POINTS: 

•CENTERBODY - 40 PERCENT FROM WEDGE L.E. EDGE TO INLET ENTRANCE. 

• COWL - HALFWAY BETWEEN INLET ENTRANCE AND THROAT. 



X/XREF 



THE COMPUTATIONAL GRID 



• GRID SIZE WAS 221 X 91. 

• GRID WAS NONUNIFORM IN X DIRECTION: 

• PACKED ON BOTH ENDS FROM THE WEDGE L.E. TO THE COWL L.E. 

• GEOMETRICALLY STRETCHED FROM THE COWL L.E. TO THE EXIT OF 
THE INLET. 

• IN Y DIRECTION. THE GRID WAS PACKED USING HYPERBOLIC TANGENT 
FUNCTION . YPLUS WAS APPROXIMATELY 1 AWAY FROM BOTH WALLS. 


• A SEPARATE GRID WAS MADE FOR THE LAUNDER AND SPALDING MODEL 
AND YPLUS OF APPROXIMATELY 30 WAS USED AWAY FROM THE WALLS. 


SCHEMATIC DIAGRAM OF P8 INLET MODEL 




BOUNDARY CONDITIONS 


NON-REFLECTIVE NO-SLIP ISOTHERMAL 



DENSITY CONTOURS FOR P8 INLET 




PRESSURE CONTOURS FOR P8 INLET 


CONTOUR LEVEL 



MACH NUMBER CONTOURS FOR P8 INLET 


CONTOUR LEVEL 
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Y/XREF 
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CONCLUSIONS 

• A COMPUTATIONAL STUDY HAS BEEN CONDUCTED TO EVALUATE 
THE PERFORMANCE OF VARIOUS TURBULENCE MODELS. 

• THE THOMAS MODEL COMPARES VERY WELL WITH THE EXPERIMENTAL 
DATA, AND IT PERFORMS BEST AMONG THE ZERO-EQUATION MODELS. 

• THE BALDWIN-LOMAX MODEL AND ITS COMBINATION WITH THOMAS 
MODEL ARE NOT ABLE TO RESOLVE THE PROBLEM OF SHOCK 
WAVE AND BOUNDARY-LAYER INTERACTION ACCURATELY. THE 
BALDWIN-LOMAX MODEL PREDICTS SEPARATION NEAR THE 
INTERACTION OF THE COWL SHOCK WITH THE WEDGE BOUNDARY 
LAYER, WHERE NONE IS KNOWN TO EXIST IN EXPERIMENTS. 

• THE CHIEN AND SPEZIALE MODEL COMPARE VERY WELL WITH THE 
EXPERIMENTAL DATA. AND PERFORMS BETTER THAN THE THOMAS 
MODEL, PARTICULARLY NEAR THE WALLS. THE LAUNDER AND 
SPALDING MODEL DOES NOT PERFORM AS GOOD AS THE CHIEN AND 
SPEZIALE MODELS. 

• AS THE CPU TIME REQUIRED FOR THE THOMAS MODEL IS FAR LESS 
THAN THE TWO-EQUATION MODELS. IT IS CONCLUDED THAT THE THOMAS 
MODEL IS BEST SUITED FOR THE PREDICTIONS OF PRESSURE 
DISTRIBUTIONS, AND THE CHIEN AND SPEZIALE MODELS ARE BEST TO 
CALCULATE FLOW QUANTITIES NEAR THE WALLS. 


PITOT PRESSURE AND TOTAL TEMPERATURE 
DISTRIBUTIONS AT THE THROAT 
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AN ALGEBRAIC TURBULENCE MODEL 
FOR TURBOMACHINERY 


by 


Rodrick V* Chima 
NASA Lewis Research Center 
Cleveland, Ohio 


OVERVIEW 


MOTIVATION - TURBINE ENDWALL HEAT TRANSFER 

. description of new model 
, RESULTS 

1. FLAT PLATE 

2. ANNULAR TURBINE CASCADE 

3. TURBINE ENDWALL HEAT TRANSFER 

4. SUPERSONIC COMPRESSOR BLADE 

» SUMMARY 


500 Rec / 1000 


EXPERIMENTAL ENDWALL STANTON NUMBER CONTOURS 
AS A FUNCTION OF 5 inlet AND Re chord 



RVC3D fROTOR VISCOUS CODE 3-D) 

BY R. V. CHIMA 
DESCRIPTION 

• EULER OR NAVTER-STOKES ANALYSIS 

FOR STEADY 3-D FLOWS IN TURBOMACHINERY 

FEATURES 

• CARTESIAN FORMULATION, ROTATION ABOUT X-AXIS 
’RECTANGULAR OR ANNULAR GEOMETRIES 

• SOLVES NAVTER-STOKES EQUATIONS 

THIN-LAYER FORMULATION, (NO STREAMWISE VISCOUS TERMS) 
RETAINS HUB-TO-TIP & BLADE-TO-BLADE VISCOUS TERMS 
BALDWIN-LOMAX OR CEBECI-SMITH TURBULENCE MODEL 
SIMPLE TIP CLEARANCE MODEL 

. NODE-CENTERED FINITE-DIFFERENCE FORMULATION 
EXPLICIT 4-STAGE RUNGE-KUTTA TIME-MARCHING SCHEME 
2ND + 4TH ORDER ARTIFICIAL VISCOSITY, EIGENVALUE SCALING 
VARIABLE Atij fe IMPLICIT RESIDUAL SMOOTHING 
HIGHLY VECTORIZED & AUTOTASKED FOR CRAY Y-MP 

• STACKED C-TYPE GRIDS 







TTIKBULENT viscosity profile 



INTERMITTENT 

TURBULENCE 

OUTER LAYER 
INNER LAYER 


rF,RF,f!T.SMTTH & B A LDWIN-T >OM A X MODELS 

TNNER LAVER: PRANDTL-VAN DR IEST - PPRMULATIQ Ij 

CEBECI-SMITH 

A i. = pl 2 \du/dy\ 

l = KyD 

D = l-exp (-y + /A + ) VAN DRIEST DAMPING 

BALDWIN-LOMAX 

Pi = P 1 * M 


OUTER LAVER: CLAUSER FORM ULATION 


CEBECI-SMITH 


Ho = KpyS'ue 


7 = 


1 + 5 


• 5 (!) 


-1 


KLEBANOFF INTERMITTENCY FUNCTION 


BALDWIN-LOMAX 

• I Vmazfmaz 

Pt> — I PI cp | wake option 

f{y) = yW\ D 


( 


BALDWIN-LOMAX MODEL ANALYSIS 

(SEE PAPER FOR DETAILS) 


1. ASSUME SUBLAYER-WALL-WAKE VELOCITY PROFILE 

2. CALCULATE BALDWIN-LOMAX FUNCTION f(y) 

MAX. OCCURS AT y mot = .6465 
INDEPENDENT OF PRESURE GRADIENT 

NO MAX. FOR INFINITELY FAVORABLE dp/dx 

3. SPURIOUS MAX. CAN OCCUR AT EDGE OF VISCOUS SUBLAYER 
MOST LIKELY AT LOW Re & FAVORABLE dp/dx 




1. VELOCITY PROFILES, Rei = 105,000 2. B-L FUNCTION /(y) 3. SPURIOUS MAXIMUM IN /(y) 


PROPOSED TURBULENCE MODEL 

TNNER LAYER (SIMILAR TO BALDWIN-LOMAX) 

Pi = M 
/ = nyD 

D = 1 - exp(-y + /A + ) 








PROPOSED TURBULENCE M QBEL 


PRESSURE GRADIENT EFFECTS 

• ACCELERATING FLOWS TEND TO RELAMIN ARIZE 

• MODELLED BY INCREASING A + IN FAVORABLE dp/ds 

• CEBECI’S EXPRESSION FOR A + USED: 


A + = 


26 

vi + rnsp* 


p + 


v dp 
pu* 3 ds 


• PRESSURE GRADIENT EVALUATED USING: 


dp~V e 

ds ~ |V,| 


• Vp 


• “EDGE VELOCITY” V t EVALUATED AT A GRID LINE FAR ENOUGH 
FROM THE WALL TO GIVE THE GENERAL FLOW DIRECTION 

• KAYS-MOFFATT EXPRESSION WAS TESTED, EFFECTS TOO STRONG 


PROPOSED TURBULENCE MODEL 

LOCAL SHEAR MODEL 

• IN STRONGLY ACCELERATING FLOWS r + DECREASES WITH y + 

• MODELLED BY REPLACING t wM WITH r(y) IN D 

D = 1 - exp(-y + /A + ) 


• ERROR IN ORIGINAL PAPER - USED p,M ONLY 

. USED BY KAYS, PATANKAR-SPALDING, OTHERS 

• ALSO USED TO AVOID PROBLEMS AT SEPARATION WHEN r waU -» 0 


PROPOSED TURBULENCE MODEL 


Ho 


7 

Cut 

Ckm 


K/y yrain 


F 

CUy(|K»,[ - \V min \) 


Jl + 5.5 

0.825 

0.55 



PROPOSED TURBU LENCE MODEL 

OUTER L AYER - FUNCTION F 
. DEFINE F = I fdy 

• INTEGRATE BY PARTS ASSUMING |wj -► 0 AS y -* 6 


F = /“yM^y 



F = 6 m u c 

. USE F DIRECTLY IN CEBECI-SMITH OUTER FORMULATION 
. ELIMINATES CONSTANT 

• DOES NOT REQUIRE KNOWLEDGE OF 6 OR u e 

. DISCOVERED INDEPENDENTLY BY D. A. JOHNSON, AIAA 92-0026 


proposed ttjr.bulencE 31QDEL 


BALDWIN -LOMAX FUNCTION f(y) 



OUTER T . AYER - LENGTH SCALE.’! 

. y IS THE CENTROID OF THE f(y) CURVE 

H fiv)dy = jl iiv)dy 

• EVALUATE USING COLE’S VELOCITY PROFILES 

n y/6 

0 .5 

.5 .55 
oo .606 

• USE EQUILIBRIUM VALUE Ckm = vh 5 = 55 


PROPOSED TURBTJ T iEN CE MODEL 

OUTER LAYER • WAKE MODEL 

^ = A/Uminj c wk y(\V maT \-\V mtn \) 

• LOWER OPTION IS A CONVENTIONAL WAKE MODEL 
. EVALUATE C wk BY EQUATING TWO OPTIONS, ASSUMING 

Vsep = -606 <5 
F scp = u e 6/2 
AV/u c « 1 


< 


• GIVES C wk = 0.825 



PROPOSED TURBULENCE MODF-T, 

3 -D IMPLEMENTATION 

• GRANVILLE BLENDING FUNCTION 

/ i,- 

M«// = Mo tanh — 

Po 

• MODEL APPLIED INDEPENDENTLY IN BLADE-TO-BLADE (77) 
AND SPANWISE (() DIRECTIONS 

• INNER LAYER - USE BULEEV LENGTH SCALE 

2 M< 

s, + S C + Js* + s 2 

• OUTER LAYER - USE ACTUAL DISTANCE ACROSS PROFILE 

t/o = Sq OR 

• BLEND 17 AND ( PROFILES VECTORALLY 

Mfuri = \Jtfq + M(< 


FLAT PLATE VELOCITY PROFILES 



COMPARISON OF FLAT PLATE VELOCITY PROFILES 
TO SPALDING’S COMPOSITE LAW OF THE WALL 







STATIC PRESSURE DISTRIBUTION 



FRACTION OF AXIAL CHORD 


COMPUTED & MEASURED PRESSURE DISTRIBUTIONS 
FOR THE ANNULAR TURBINE CASCADE 



EXIT P 0 LOSS COEFFICIENT 



COMPUTED & MEASURED LOSS COEFFICIENT PROFILES 
FOR THE ANNULAR TURBINE CASCADE 






turbine cascade 



// 
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SUMMARY 


. SPURIOUS MAXIMUM IN B-L FUNCTION f(y) CAN ^COR- 

RECT TURBULENT LENGTH SCALE fc ERRATIC St OR Cj PATTERNS 


- MOST LIKELY AT LOW Re AND FAVORABLE dp/ds 

• NEW TURBULENCE MODEL PROPOSED 

- INTEGRAL RELATIONS FOR 6*u e AND <5 USED WITH C-S MODEL 

- EFFECTS OF dp/ds MODELED 

- WAKE MODEL PROPOSED 

• FLAT PLATE 

- B-L & NEW MODEL AGREE WITH LAW OF THE WALL 

- LOCAL SHEAR MOD. DOES NOT AGREE WITH LAW OF THE WALL 

• ANNULAR TURBINE 

- GOOD AGREEMENT WITH EXPT. PRESSURE DISTRIBUTION 

- WAKE MIXING UNDER-PREDICTED 


* TURBINE ENDWALL HEAT TRANSFER 

- VARIATIONS IN ENDWALL St WITH Re PREDICTED WELL 

- EFFECTS OF dp/ds IMPORTANT 

• TRANSONIC FAN 

- SHEAR LAYER FROM BOW SHOCK ACTS LIKE VISCOUS LAYER 

- NEW MODEL OVERPREDICTS L.E. p t 

- B-L MODEL PREDICTS REASONABLE L.E. 


Lewis Research Center 

National Aeronautics and 
Space Administration 


PROPULSION SYSTEMS DIVISION 
COMBUSTION TECHNOLOGY BRANCH 



LOW EMISSIONS COMBUSTORS 


J. M. DEUR 

SVERDRUP TECHNOLOGY, INC. 
LEWIS RESEARCH CENTER GROUP 
BROOK PARK, OHIO 


Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


APPLIED ANALYTICAL COMBUSTION/EMISSIONS RESEARCH 

• ANALYZE LOW EMISSIONS COMBUSTORS TO AID IN-HOUSE EXPERIMENTS 
AND CONTRACTOR COMBUSTOR DEVELOPMENT PROGRAMS. 

• PRESENT: 

• UTILIZE EXISTING CODES , PRINCIPALLY KIVA-ll. 

• IMPROVE NUMERICS AND PHYSICAL MODELS (E.G., PDF COMBUSTION- 
TURBULENCE INTERACTION) ON LIMITED BASIS TO SATISFY 
CRITICAL NEEDS . 

• FUTURE: 

• ADOPT 3-D ALLSPD CODE WHEN AVAILABLE. 




Lewis Research Cen.er PROPULSION SYSTEMS DIVISION 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


KIVA-II FEATURES 

MULTI-DIMENSIONAL TIME ACCURATE FINITE DIFFERENCE CODE. 
COMPRESSIBLE FLOWS. 

, k-e TURBULENCE MODEL WITH WALL FUNCTIONS OR SUB-GRID SCALE 
TURBULENCE MODEL. 

. AD iabATIC OR CONSTANT TEMPERATURE WALL BOUNDARIES. 

. ARBITRARY MESH. 


Lewis Research Center 

National Aeronautics and 
Space Administration 


PROPULSION SYSTEMS DIVISION 
COMBUSTION TECHNOLOGY BRANCH 


APPLIED ANALYTICAL COMBUSTION/EMISSIONS RESEARCH TEAM 


University of Florida (I) 
Prof Jerry Mickiow 
(KIVA-II) 


Michigan Tech University 
Prof. Jason Yang 
(KIVA-II) 


Los Alamos National Laboratory 
Dr. Mike Cline 
(KIVA-II) 


NASA Lewis Research Center 
Drs. John Deur/Krish Kundu 
(KIVA-U/LSENS) 


Carnegie Mellon University 
Profs. Tom Shih/Juan Ramos 
(LeRC-3D) 


University of Florida (II) 
Prof. Nick Winowich 
(KIVA-II) 


University of South Florida 
Prof. Ben Ying 

(KIVA-II) 

Carnegie Mellon University and University of South 
Florida leaving program at end of FY93 
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RICH BURN - QUICK MIX - LEAN BURN (RQL) FLAME TUBE 






Lewis Research Center 

^r.=^ and COMBUSTION TECHNOLOGY BRANCH 


PROPULSION SYSTEMS DIVISION 


RQL FLAME TUBE WALL TEMPERATURE COMPARISON 


iir. 



SLANTED SLOT MIXER (W/0 SWIRL) TWO HOLE MIXER (W/0 SWIRL) 





ONE HOLE MIXER (W/0 SWIRL) 



M: 



ONE HOLE MIXER (W/ SWIRL) 



Lewis Research Center 

PROPULSION SYSTEMS DIVISION 

K 

National Aeronautics and 
Space Administration 

COMBUSTION TECHNOLOGY BRANCH 
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RQL FLAME TUBE TEMPERATURE COMPARISON 






Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 



RQL FLAME TUBE NO x EMISSION INDEX COMPARISON 



ONE HOLE MIXER (W/O SWIRL) ONE HOLE MIXER (W/ SWIRL) 



MIN NOTE : CALCULATION CONSIDERS THERMAL NO x ONLY. MAX 


Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


RQL MIXER CONFIGURATION PARAMETRIC STUDY NO x COMPARISON 


SLANTED SLOT MIXER W/0 SWIRL TOO 

TWO HOLE MIXER W/0 SWIRL 1.25 

ONE HOLE MIXER W/0 SWIRL 1.39 

ONE HOLE MIXER W/ SWIRL 0. 75 

SLANTED SLOT MIXER W/O SWIRL (EXPERIMENTAL) TOO 


NOTE: 

• VALUES ARE SPATIAL AVERAGES TAKEN AT SAMPLING LOCATION B. 

* VALUES ARE NORMALIZED BY SLANTED SLOT MIXER EXPERIMENTAL 
READING. 
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THE EFFECTS OF TURBULENCE MODELING ON THE 
NUMERICAL SIMULATION OF CONFINED SWIRLING FLOWS 


G. J. MICKLOW AND M. R. HARPER 
UNIVERSITY OF FLORIDA 
GAINESVILLE, FLORIDA 

J. M. DEUR 

SVERDRUP TECHNOLOGY, INC. 
BROOK PARK, OHIO 


AIAA/SAE/ASME/ASEE 29TH JOINT PROPULSION CONFERENCE 
MONTEREY, CALIFORNIA 
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K - £ TURBULENCE MODEL 


DK _ d [ ■ u, 3K \ 3uj / 3uj 3u| 
Dt 3xj [ 3xj j * 3xj I 9xj 3xj 


C K 


C M = 0.09 




Lewis Research Center 

PROPULSION SYSTEMS DIVISION 

i 

National Aeronautics and 1 

Space Administration 

COMBUSTION TECHNOLOGY BRANCH 



K - £ TURBULENCE MODEL 


Dt 3xj l o t - 9xj 


K 3X; 1 0X; 9x; 


C 1 =1.44 
C 2 -192 
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C„ MODIFICATIONS TO K - £ TURBULENCE MODEL 

C U = ^ 

“ k" 2 

A = 0.0083 


C u = 0.09 
V wall 
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C 2 MODIFICATIONS TO K - £ TURBULENCE MODEL 


:C 2 (1 - C3 Ri) 
= 1.92 
* 0.2 


(0.1 <c 2 <2.4) 


C 2 =1.92 
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UCI AXISYMMETRIC CAN COMBUSTOR 
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COMPARISONS TO UCI AXISYMMETRIC CAN COMBUSTOR DATA (z=24 cm) 


* oase 


..ocoi>o« -* c>" 


VELOCITY CORRELATION 


*)ad>o! LOfutiCf' ' «. err 


TEMPERATURE 
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CONCLUSIONS 

• INLET VELOCITY IS BEING MODIFIED TO MATCH FIRST INLET STATION. 

. PDF COMBUSTION-TURBULENCE MODEL OF HSU, ET AL. IS BEING ADDED. 

• MUCH MORE WORK IS REQUIRED. 
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DEVELOPMENT OF A 
RELIABLE ALGEBRAIC TURBULENCE 
MODEL GIVING ENGINEERING ACCURACY 
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by 


B.P. Leonard and J.E. Drummond 
The University of Akron 
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I. CHOICES FOR TURBULENCE MODELLING 

A. ZERO-EQUATION (ALGEBRAIC MODELS) 

B. MULTIPLE-EQUATION MODELS 


II. THE NEED FOR A RELIABLE ALGEBRAIC MODEL 

A. RELATIVE SIMPLICITY 

B. ENGINEERING ACCURACY 

C. COST-EFFECTIVE APPLICATION 

III. DEVELOPMENT OF THE MODIFIED MIXING LENGTH (MML) 
MODEL 

IV. COMPARISON OF MODELS 

A. MML 

B. BALDWIN-LOMAX 

C. TWO-EQUATION (k-e) MODEL 


V. WHERE DO WE GO FROM HERE? 
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EFFECTIVE VISCOSITY: fi eff = [i + /*, (1) 

TURBULENT VISCOSITY: n, = P ? 2 M (2) 


MIXING LENGTH: 


(CA . 
— My 

i - 

(i - 

1 - exp ( y )1 for y + < C, 

(3) 

\ r V 


\ Cj/ . 

L \ A + / J 



i 

- ■(§)’• 

for y * > 

(4) 


where: y + = y/y*; y* = p/\[p I 
SHEAR STRESS “FILTER”: 

|r,| = 0.1 |r { _ 2 | + 0.2 | t._, | + 0.4 |rj + 0.2 \r M \ + 0.1 |r M | (5) 




LeRC 

CTM 9/93 
WORKSHOP 


REFINING MML - CAPPING LENGTH SCALE AS 

A FUNCTION OF POSITION 


■ 


OUTER REGION LENGTH SCALE (CAPPING LENGTH): 

£ = £ = K - V * 

* * cap * q y > 


CAPPING LENGTH BASED ON BOUNDARY LAYER THICKNESS: 

* cap = B.S 


CAPPING LENGTH BASED ON LOCAL SHEAR STRESS: 


C AP = (3.31 • 10' 6 )B | C f | " 3 - 5 y * 
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VARIATION OF MIXING LENGTH 
IN ATTACHED BOUNDARY LAYER 
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SHEAR LENGTH SCALE, y*, 
BASED ON |r w | AND EQUATION (5) 
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THEORETICAL BASIS OF 
BALDWIN-LOMAX MODEL 


INNER REGION TURBULENT VISCOSITY: 


WHERE 


\i t = p i 1 | co 


i = K y[l - exp (~y*IA + )] 



OUTER REGION TURBULENT VISCOSITY ASSUMES THE SMALLER 
VALUE OF: 

= P K C cp F k (y) y max F max 


p, = 0.25 p K C U 2 m y m JF n 


WHERE 


F(y) = y | co 


c k y\ 6 r 


0.00050 
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FLAT PLATE GEOMETRY AND GRID 
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UNIFORM 
MACH 0.2 
FLOW 


/ 


DEVELOPING BOUNDARY LAYER 
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MML FLAT PLATE SKIN FRICTION 
CALCULATIONS FOR VARYING C { /8 + 


O WEIGHARDT DATA 

MMIT (Cl/6*- .3) 

- MMLT (Cl/6** .4) 

• MMLT (Cl/6*- .5) 

■ - MMLT (Cl/6*- .6) 
- -MMLT (Cl/6*- 7) 
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MML FLAT-PLATE VELOCITY DEFECT 
PROFILES FOR VARYING C,/6 + 

13 
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SKIN FRICTION CALCULATIONS FOR 
VARIOUS TURBULENCE MODELS 


0.0035 ro/jv 
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0.0025 

'f 

0.002 

0.0015 

0.001 

0.0005 

0 


0 ' 

o WEiGHARDT DATA 

MML.T (C1/6*« .7) 

■ * • * BAIDWIN-LOMAX 
K-EPSILON 


6000 p^ 8000 
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VELOCITY DEFECT PROFILES FOR 
VARIOUS TURBULENCE MODELS 
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TURBULENT VISCOSITY PROFILES 
FOR VARIOUS MODELS 
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CONVERGENCE CHARACTERISTICS FOR 
FLAT PLATE CALCULATIONS 


MODEL 


ITERATIONS 


CRAY Y/MP 
CPU TIME (s) 


MMLT 

BaJdwin-Lomax 
Chien k-e 


M .. 

rtuio . 



















GEOMETRY FOR BACKWARD 
STEP CALCULATIONS 
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GRID FOR BACKWARD STEP CALCULATIONS 
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INCORPORATING THOMAS MODEL 
FOR FREE SHEAR LAYER 



LENGTH SCALE FOR UNBOUNDED REGION FOR THOMAS MODEL: 


_ f o [Max(| Uj |) - Min(j Uj |)J 


TURBULENT VISCOSITY IN TRANSITION REGION BETWEEN 
BOUNDARY LAYER AND UNBOUNDED REGION: 


P tr = 


^MML (^4 y ) + ^Th (y ^"3 ) 


c - c 

^4 '"'3 
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WALL BOUNDARY AND FREE SHEAR LAYERS 
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BACKWARD-FACING STEP SKIN FRICTION 
CALCULATIONS FOR VARIOUS MODELS 
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REATTACHMENT POSITION FOR 
VARIOUS TURBULENCE MODELS 



CASE 


Driver-Seegmiller Data 

MMLT (C1/6+ = .7, 

C3 = 10£+cap.C4 = 20^*cap) 

Thomas 

Baldwin-Lomax 
Chien k-E 


REATTACHMENT 

posmoN 

(STEP HEIGHTS, H d ) 


6.25 

7.42 
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CONVERGENCE CHARACTERISTICS FOR 
VARIOUS TURBULENCE MODELS 
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WHERE DO WE GO FROM HERE? 



A. APPLICATION TO MORE COMPLEX FLOWS 

• NOZZLE FLOWS 

• RECIRCULATING INTERNAL FLOWS 

• EXTERNAL FLOWS 


B. COMPARISON WITH OTHER TURBULENT CALCULATIONS 


C. ESTABLISH REASONABLE ALGEBRAIC MODEL AS AN OPTION 
TO GIVE RESULTS WITH ENGINEERING ACCURACY FOR A 
VARIETY OF APPLICATIONS 
















Application of Algebraic and 
Two-Equation Turbulence Models 
to HSR Nozzle Flow Calculations 


N. J. Georgiadis 
and 

J. R. DeBonis 

Nozzle Technology Branch 
Propulsion Systems Division 


September 16, 1993 


FULL NA VIER-ST OKES ANALYSES 
OF NOZZLES (WITH PARC CODE) 


OVERVIEW OF PARC CODE AND AVAILABLE TURBULENCE 
MODELS 


VALIDATION TEST CASES: 

1. AXISYMMETRIC PLUG NOZZLE 

2. EJECTOR NOZZLE 

HIGH-SPEED RESEARCH (HSR) NOZZLES WITH 
POTENTIAL FOR NOISE REDUCTION: 

1. NASA/GE 2DCD MIXER/EJECTOR NOZZLE 

2. PRATT & WHITNEY 2D MIXER/EJECTOR NOZZLE 



OVERVIEW OF PARC 


• 3D AND 2D/AXISYMMETRIC VERSIONS 

• NAVIER-STOKES AND EULER MODES 

• CENTRAL DIFFERENCE DISCRETIZATION 

• BEAM AND WARMING ALGORITHM 

• SEVERAL TURBULENCE MODEL OPTIONS 

• CAPABILITY TO HANDLE MULTIPLE GRID BLOCKS 

(NONCONTIGUOUS INTERFACING) 

• GENERALIZED BOUNDARY CONDITIONS 


Full Navier-Stokes 


at 


' p " 


P U i 

E 

JT1 

II 

i 


Equations in PARC 


= 1 3g 3 

9Xj Re dxj 


p Uj 


0 

>J + P8 u 



+ P)Uj 


. u * T jk - Qj . 




5Uj 

5X: 3Xj 

V J 

-k, aT 


2 9u k ~ 


^ = 


(Y-l)Pr ^ 


M'eff = fl + Hi 


K 


c-eff 


= K c + 




Thomas Model 


ii t = p^ 2 H 


Wall bounded flows: 


f 

£ = Ky 1 



\ 


V 2 


£ o Max^UjJj - Min(|uJ] 
co c 


Free shear layers: 


Baldwin-Lomax Model 


Pt = 


[ (P’t ) inner ’ y — Yc 

[(^tWr . y ^ y 0 


(Pt)outer " p[^C CP F waIce ] Fj 


kJeb 


^wake (ymax^max ’ ^wkymax U dif ^ Fjnax ) 


F(y) = y |co| 


1 - e 


-y 

A* 


F kleb(y) = 


1 + 5.5 


r C v^ 6 

^kleb Y 
V y max J 


-i-I 


Chien k-8 Model 


= C.f, p k 2 / £ 


D(pk) 

Dt 


_a_ 

3x s 


( p-t 1 

3k 

p+— 

.V °k J 

9xj 


+ n-pe-D 


D(pe) 

Dt 


3 

Pil 

ii" 

3x ; 

p+- £ - L 

A **£ > 

3Xj 


+ C a f,ni - C £ 2 f ; p|- - 2 p^-e' 


(-0.5y + 


D = 2|ik/y 


= 1 - e 


(-.0115 y*) 


f,= 1.0 


*2 


1.0 - 0.22 e 


-(R,/6) J 


Re 



two - fottation model of speziale 


• SOLVED AS A k-e MODEL IN PARC 

• ORIGINALLY DEVELOPED AS A k-T MODEL 
IN NASA CR-182068 

• UPWIND DIFFERENCE DISCRETIZATION 
(MAIN FLOW EQUATIONS IN PARC USE 
CENTRAL DIFFERENCING) 



BENCHMARK VALIDATION CASES 


1. Mach 0.2 flat plate 

2. Mach 0.128 backward-facing step 


SHEAR STRESS PREDICTIONS WITH THE 
DIFFERENT TURBULENCE MODELS IN PARC 





SHEAR STRESS AFTER ™ESTEP WITHTHE 
DIFFERENT TURBULENCE MODELS IN PARC 



LANGLEY SINGLE FLOW PLUG NOZZLE 


• VENTED AND NON-VENTED PLUGS 

• 15° PLUG HALF ANGLE 

• HEAVILY INSTRUMENTED TO MEASURE: 

i p L xjG SURFACE TEMPERATURES, PRESSURES, SHEAR STRESS 

2. JET PLUME QUANTITIES ( INCLUDING LDV & FLOW VISUALIZATION) 

3. FLOWFIELD ACOUSTICS 




C L 


PRIMARY FLOW NOZZLE 



GEOMETRY FOR GRID GENERATION 
AND PARC2D CALCULATIONS 



VELOCITY PROFILES FOR FINE GRID SOLUTIONS 


X = 25 in. 



X = 50 in. X = 75 in. 




0 500 1000 1500 2000 2500 3000 

V (ft/s) 


500 1000 1500 2000 2500 3000 

V (It/S) 


500 1000 1500 2000 2500 3000 
V (ft/s) 




TOTAL 

FOR 

TEMPERATURE PROFILES 
FINE GRID SOLUTIONS 


X = 25 in. 

1 f. -j. 

X = 50 in. 

X = 75 in. 
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Thomas 
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K-Epsilon 


Total Temp. (°R) 


Total Temp. (°R ) 


Total Temp. ( R) 


SHOCK FUNCTION 
(BASED ON PRESSURE GRADIENT) 
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/Thomas 
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K-Epsilon 



STATIC PRESSURE DISTRIBUTION ALONG PLUG 



2D FTECTOR NOZZLE TEST CA SE 


»p'PPF < sFNTA r nVE OF MIXER-EJECTOR NOZZLES 
THAT^AEJE^ BEING CONSIDERED FOR SUPERSONIC 
TRANSPORT APPLICATION 

FLOW DOMINATED BY TURBULENT MIXING OF A 
HIGH ENERGY STREAM WITH SECONDARY AIR 





2D EJECTOR NOZZLE OF 
GILBERT AND HILL 













TURBULENT VISCOSITY CONTOURS 
FOR 2D EJECTOR NOZZLE 
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AXISYMMETRIC INVERTED VELOCITY 
PROFILE NOZZLE OF VON GLAHN, 
GOODYKOONTZ, AND WASSERBAUER 



AXISYMMETRIC UNDEREXPANDED NOZZLE 
OF HELSTLEY AND CROSSWY 



Mj e , = 1.56 

M 00 = 0.60 

NSPR = 4.5 



Mixer/Ejector Nozzles 

• Entrain large amounts of secondary flow 

• Rapidly mix two flows together to lower 
jet velocity 

• Lower jet velocity results in lower noise 

• Maintain high thrust due to large mass 
augmentation 


F = mv 


Typical Mixer/Ejector Nozzle 
Chute Geometry 



Secondary Flow 


NASA/GE 2DCD 
Mixer/Ejector Nozzle 





Shroud Static Pressures 


MAR = 1.20 

Thomas Turbulence Model 



X/L 




Shroud Static Pressures 
MAR = 1.20 










Flow Angle Contours 
Nozzle Exit Plane 



PARC3D K- e PARC3D Thomas 


Anatomy of the HSR P&W 2-D Mixer Ejector Nozzle 


Inventory: 

• 2 Mixer Designs 

• 3 Shroud Lengths 

• 3 Shroud/Sidewall 
Acoustic 
Treatments 

• Sidewalls with 
Windows for Flow 
Visualization 


Mixer 



Sidewalls 





P&W MIXER/EJECTOR NOZZLE 
Mach Number Contours 
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Peak Side 
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2nd GENERATION MIXER EJECTOR ANALYSIS 


Exit Plane Total Temperature 


UTC/NASTAR 



CHOl/PARC 




Conclusions 

• Algebraic models are inadequate for 
predicting mixer nozzle flowfields 

• Two-equation model available in PARC3D 
give no improvement over algebraic models 

• A reliable two-equation model is required 
in PARC3D for complex nozzle flows 





Turbulence Modeling in Aircraft Icing 

N94-18555 


Mark G. Potapczuk 

Icing & Cryogenic Technology Branch 


Workshop on Computational Turbulence Modeling 
NASA ICOMP 
September 15, 1993 


INTRODUCTION 


The Icing and Cryogenic Technology Branch develops computational tools 
which predict ice growth on aircraft surfaces and uses existing CFD 
technology to evaluate the aerodynamic changes associated with such 
accretions. 


Surface roughness, transition location, and laminar, transition, or turbulent 
convective heat transfer all influence the ice growth process on aircraft 
surfaces. 


Turbulence modeling is a critical element within the computational tools 
used both for ice shape prediction and for performance degradation 
evaluation. 



CURRENT CODE DEVELOPMENT 


2D CODES 


LHVICE - POTENTIArFU^QNTE^BffiBSBNpAlX LAYEF 


LEWICE/IBL - POTENTIAL FLOW / INTERACTIVE BOUNDARY LAYER 


LEWICE/NS - NAVIER-STOKES; STRUCTURED GRID 
LEWICE/UNS - NAVIER-STOKES, UNSTRUCTUREP;GRID| 


3D CODES 

. LEWICE3D - PANEL CODE / INTEGRAL BOUNDARY LAYER 
• LEWICE3DGR - ANY GRID BASED FLOW SOLUTION 


ICE ACCRETION MODELING 

CURRENT MODEL USED FOR ICE GROWTH 


. MASS AND ENERGY BALANCE IN CONTROL VOLUMES ALONG THE 
SURFACE 


. CONVECTIVE HEAT TRANSFER IS MAJOR FACTOR IN ENERGY 
BALANCE 

. INTEGRAL BOUNDARY LAYER FORMULATION USED TO DETERMINE 
LAMINAR AND TURBULENT HEAT TRANSFER COEFFICIENTS 


. SURFACE ROUGHNESS MODELED 
ACTUAL ICE ROUGHNESS VARIES 
BOUNDARY LAYER THICKNESS 


AS SAND-GRAIN ROUGHNESS; 
FROM SMALLER TO LARGER THAN 


ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 


SKIN FRICTION COEFFICIENT 


^ = 0.1681 


In 


/864.O0 


- + 2.568 




i-2 


V 


) 


WHERE 


e ,w = 


i0.8 


0.0156 


v 4 


ii 


tr 


V™ 6 ds 

e 


+ w 


ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 
LAMINAR 



- 1/2 


TURBULENT 


h t (s) = StpV e c p 


cJ 2 

Pr t + \c/2(l/St k ) 


P V A 


ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 


ROUGHNESS STANTON NUMBER 



Vk 

1.16 (H“) 


- 0.2 


AND 



ICE ACCRETION MODELING 

ICE ROUGHNESS CHARACTERIZATION 
SAND-GRAIN ROUGHNESS 


ACTUAL ICE ROUGHNESS 



ICE ACCRETION MODELING 


PLANS 


• EXPERIMENTS TO CHARACTERIZE ICE ROUGHNESS GEOMETRIES AT A 
VARIETY OF ICING CONDITIONS 


• EXPERIMENTS TO CHARACTERIZE VELOCITY FIELD OVER REAL AND 
ARTIFICIAL ICE ROUGHNESS GEOMETRIES 


• EXPERIMENTS TO MEASURE HEAT TRANSFER OVER REAL AND 
ARTIFICIAL ICE ROUGHNESS GEOMETRIES 


. DEVELOPMENT OF MODIFIED COMPUTATIONAL MODEL BASED ON 
THESE EXPERIMENTS 


ICED AIRFOIL AERODYNAMICS 


NACA 0012 ICING CONDITIONS 


a = 4° V = 130 mph 

d = 20 ja.m LWC = 2.1 g/m 3 
T = 1 8° F 


SHEAR LAYER 



) 



ICED AIRFOIL AERODYNAMICS 

baldwin-lomax turbulence model 



/ = mixing length 



LL ~F Y 

f* max 


max max 


F (y) = y|co| ( 1 - exp ( (-/) / A) ) 


ICED AIRFOIL AERODYNAMICS 

BALDWIN-LOMAX TURBULENCE MODEL 


NORMAL B.L. F PROFILE RECIRCULATION REGION 

F PROFILE 




ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 

H , = P^ 2 ® 

WHERE 


y+ty 

* 

y 


< ^1 


*00 




ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 



ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


THE CEBECI-CHANG ROUGHNESS MODEL IS ADDED TO THE 
TURBULENCE MODEL 



WHERE, 

A y + = (40 (w T /v) and k* = k s (u % /v) 


ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 



AOA 


ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 



ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


STRUCTURED GRID FOR ARTIFICIAL ICE SHAPE 



0.00 0.05 0.10 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


UNSTRUCTURED GRID FOR ARTIFICIAL ICE SHAPE 



ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 
STRUCTURED GRID MACH NUMBER CONTOURS 



Q 

R* 

GRID 



ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


UNSTRUCTURED GRID MACH NUMBER CONTOURS 


in 



CONCLUDING REMARKS 


. TURBULENCE MODELING PLAYS A ROLE IN ICE GROWTH PREDICTION 
AND IN PERFORMANCE EVALUATION 


. NEW MODELING IS REQUIRED FOR THE LARGE ROUGHNESS 
ELEMENTS OF A TYPICAL ICE ACCRETION 


. AN EXPERIMENTAL PROGRAM IS CURRENTLY UNDERWAY TO 
DEVELOP A DATABASE FOR CREATION OF SUCH A MODEL 


. AN ALTERNATE ALGEBRAIC TURBULENCE MODEL HAS BEEN USED TO 
EVALUATE PERFORMANCE DEGRADATION DUE TO ICING 

. THE MML MODEL HAS BEEN USED IN AN UNSTRUCTURED GRID 
NAVIER-STOKES CODE TO CALCULATE FLOW OVER AN ARTIFICIAL ICE 
SHAPE 



Applied RNG Algebraic Turbulence Model for 
Three-Dimensional Turbomachinery Flows 


K. R. Kirtley 


Sverdrup Technology, Inc, LeRC Group 

and 

Cambridge Hydrodynamics, Inc. 


What is RNG ? 

Start — ►Navier-Stokes 


RNG 

1) Decompose velocity Into low and 
high wave number components 

2) Use perturbation theory to eliminate 

high wave number bands then renormalize 
spectrum and repeat to infrared cutoff 

3) Correlations disappear through mode 
elimination procedure 

4) Closure is automatic 

5) Evaluate coefficients from high Re 
limit (fixed point) of perturbation 
expansion 


Reynold’s Averaging 

1) Decompose velocity into mean 
and fluctuating components 

2) Ensemble (time) average over 
entire spectrum 

3) Correlations arise 

4) Model correlations (Boussinesq) 

5) Evaluate coefficients from generic 
flow data, e.g., wake, jet, etc. 

Problem Dependent 


Self Consistent 

Result — ►effective viscosity 


RNG-Based Algebraic Model 


V — Vo 


1 + H 





1 

3 


e = P = v S 


where 


S = 



duj 

dxi 




dui 

dxj 


A / = Infrared cutoff = top of the Inertial range 
resolvable < A y< modelled 


a 

C c 


= .1186 from -5/3 spectral law 

400 

=jtfUom a and the ultraviolet dissipation range cutoff 


How do you go from Fourier space to physical space ? 

why 27 r of course. 

Ti ^ 

Integral scale, Lj- — — 


But we need a mixing length, therefore 

2. _ k 

Assume E = C jf€ 3 A 3 


= 1.6075 


— ^ 2 


Integration from a, to infinity gives k = -Ck^j 3 6 

k 2 

Assume in high Re limit v = ^ 2 gh z= Cp — 

Then L/ = 4 Ay 1 


Giving 


V = Vq 


1 + H ( —z L\ — C c 


Vr 


Attributes of the RNG-Based Model 


1) The Heaviside function mimics: 


a) Near wall damping 

v — v 0 




= 9.2 


K 


b) Intermittency 

C c v* 

v = Vq when e < 

e o In the outer flow f 

c) Transition 


See Next Slide 


2) Energy-Based ^Non-equilibrium effects can be included 
through clever manipulation of the dissipation rate 


Cn 

© 


* u-v 




r P 


Cubic vs. Quartic 


y yf[~l + H 0~$- (<)j V*- vy* 1 - (rtf- vC 


Good: Analytic Solution 

Bad: Multivalued 


Single Valued 
Highly Non-linear 



S L. 4 / 1 


4 /,. 2 
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High Re limit 


4k / 


RNG model 


T LfV^o 2 




Couette Flow 


Friction Factor 

Near Wall Behavior vs. 

Reynolds Number 




Turbulence Length Scales 


In boundary layers 


C.4- 


/ K7l 

Lf = Cy 5 tanh ^ 


K and Q~ 0 0845 from RNG 

In wake region use Ra| & Lakshmlnarayana correlation for 
cascade wakes L{ = min ( KS) C w b) 

0.58 


6 = <5 0 + c CJ 1.35 — 1- 0. 

c 


(; +002 ) 


Cd = coefficient of drag ( assume 0.015) 

C = local chord 

< 5 0 = ave. of s.s and p.s. trailing edge boundary layer thickness 
C w = 0.169 from wake behind a circular cylinder 




Transverse Velocity 






















Conclusions 


# RNG based algebraic model is a viable alternative 
to traditional algebraic models for complex flows, 

# RNG model mimics transitions and near-wall 
..damping through the Heaviside function 


# Good boundary layer profiles, loss profiles, pressure 
loadings can be achieved for viscous dominated flow 

# Stanton number predictions are not yet adequate 
when using the "plain vanilla" form of the model 

# Work is needed to include non-equilibrium effects 
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Introduction 

Diffusing S -ducts common on aircraft: 

727, L-1011, F-16, F-18 

S-duct fabricated and tested 

CFD: 3-D full Navier-St okes, two turbulence models 

Research extends previous PNS and FNS studies 

Overall features qualitatively correct: total pressure 
velocity vectors, exit vortices 

Boundary layer separation details generally 
not correct 

Turbulence models and/or grid resolution usually 
cited for lack of agreement 


Code Selected: PARC3D 

Originally ARC3D - Pulliam and Sfeger 
Modified for internal flows - Cooper 

Full 3-D Navier-Stokes equations in Reynolds (mass) 
average form 

Beam-Warming approximate factorization 
Multi-block for computer efficiency 

Bow Reynolds number K - e turbulence model 
of Speziale employed - Nichols 

Algebraic Baldwin-Lomax turbulence model 




Boundary Conditions 


No slip at wall 

Total pressure and temperature specified at inlet 
Static pressure specified at exit 
S ymmetry about x - z plane 
J< = 0 on wa ll 


2H 




_dy 

wall_ 


Jt, e : zeroth -order ex trapolat ion at inlet, 
outlet, and centerline” 


Nomenclature 


C n P 0 P wall 

0 " p . p 

P 0 c i \vall 


P ~ P wall 
P 0 d P wall 


+ 

u = 


u 


U‘ 



T 

*_\ 

w 

\ 

p 

\ 

w 


*w 

ip , u f 
2 cl cl 




luc^resFconaraons am 


reoi 


Re = 2.6 x 10° 

M J= 0.6 

D 1’ D 2 = 8-04, 9 - 90 in 

4^ = 1.52 


e 


max 


12 =30 


measurement planes S/D j -0.5, 5.73 
upstream and downstream pipes =3.75Dj[ 


Grid for S-duct 


1 


exit I 

1 inlet ^1111 



1 o Grid in 3 blocks 

32x71x53 
69 x 71 x 53 
32 x 71x53 


1 H Grid in center 

129 x 11 x 15 




The geometry of the diffusing S-duct 

* 



Total pressure contours at s/D* = -0.5 
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at s/Di= -0.5 



Boundary layer wall coordinate plots 

at s/Dj = 5.73 









Streamlines along S-duct centerline 



Conclusions 

Computed flow fields agree reasonably well with 
experimental flow field 

K - e turbulence model better predicts pressure field 
than algebraic model 

Both models underpredict length, angular extent 
and axial location of boundary layer separation 

Possible causes include: inappropriate artificial/ 

computed viscosity, especially in separated region 

Improvements are needed to account for strong 
secondary flows withseparation 



Proteus Experience with the Modified MML 
Turbulence Model 


Julianne Conley 
NASA Lewis Research Center 


The Workshop on Computational Turbulence Modeling 
NASA Lewis Research Center 
September 15-16, 1993 



Overview 


• Background 

• Modification of MML Turbulence Model 

• Test Cases 

• Concluding Remarks 


Background 

• Based on original MML model of Potapczuk (1989) 

• Simple model for flows where Baldwin Lomax model falls short 

• Based on Prandtl's mixing length theory: 

M, = p/ 2 |w| 



J.O 





Modification of the MML Turbulence Model 

Step 1: Get good estimate of x w 

Step 2: Evaluate and modify MML for zero pressure gradient 
boundary layer flows 

Step 3: Modify model for an adverse pressure gradient flow 
Step 4: Combine all features into one general model (MMLPG) 





Step 1: Get good estimate of 


3w, 

x - “ ^ 


Global approach - use momentum equation and 2 interior grid points 

3x 

3y, 


= dp 
5x1 



X 


x = ay 2 + by + c 


a = 


^ l -^ 2 -b{y l -y 2 ) 


y\-y\ 


»=% 


c = ^ (Xj + x 2 — a O'} + y\) -b(y l +y 2 )) 


=> X = c 

w 


,du . 

^ ^ total 


Use local average of x w to avoid large values of y* near separated regions: 
lx. I = 0.1|x._ 2 | + 0.2|x._ j| + 0.4|x.j +0.2|x j . + 1 | + 0.1jx. +2 | 




n/(n - °°n) 


Step 2: Evaluate and modify MML for zero pressure gradient 
boundary layer flows 


Outer length scale: 


C, * 
l = , 


y*>Cx 


Doesn't allow boundary layer thickness to grow properly. 


Baldwin and Lomax 



Re x =7 x 10 6 

Re x =10x10 6 

Re x =13x10 6 

Klebanoff-Diehl 

Schultz-Grunow 


Re x =7 x 10 6 
Re x =10 x 10 6 
Re x =13x10 6 
Klebanoff-Diehl 
Schultz-Grunow 



=L ft 

^ 200 
3- L 


Turbulent Visnnsitv for Flow over a Flat PI 


Re, = 7x1 0 6 
Re x = 10x10 6 
Re, = 13x10 6 


0.5 1.0 1.5 2.0 

y/8 










Step 3: Modifications for adverse pressure gradient flow 


Benchmark test cases : Equilibrium boundary layer flows 

(P. Bradshaw, 1966) 




R = — ^ = constant 
F % dx 


a = 0 ,- 0 . 15 ,- 0.255 


P = 0,1,5 





(u_ - u)/u 


Comhininn with Blending Function: 


i + = K c\— 1 1 - ( I - ~ ) 1 - c 


kC >c7 ’ 


C, = 0 .4 C 4 , 
c 3 = /c« (p) 


C 2 = 5C 3 k 


C 4 — /c n ( P , c f ) 


C 4 - C 5 + C $ c f 


Strength 

of 

pressure 

gradient 

P 

c 3 

C 5 

c 6 

zero 

0 

1.00 

23,300 

-7.75x1 0 6 

mild 

1 

1.25 

30,100 

-1.16X10 7 

strong 

5 

1.53 

33,800 

-2.09x1 0 7 











Step 4: Combine to get one general model, MMLPG 


C 3 = 1.0 


C 5 = 23, 300. 


C 6 = -7.75 x 10 6 


C 3 = 1.0 + 0.307(3- 0.0391 (3 2 
C 5 = 23, 200 + 8560(3 - 1230(3 Z 
C = -7.75 x 10 6 - 4.51 x 10 6 (3 + 386, 000(3 2 


C, = 1.52 


C 5 = 33, 900 
C, = -20, 900 


Velnritv Defect Profiles - Mild Adverse Pressure Gradient 


MMLPG 



- X = 1.2 m 
x - 1.7 m 
x«2.1 in 

f > Bradshaw Dal a 





0 2 ii a n.i; 0 B 


Baldwin and Lomax 


— x «= 2.1 in 
o Bradshaw Data 



Static Pressure Distribution - Weak Shock Case, R=0.82 

Top Wall 


j 

MMLPG 

1 

D4«wm-uymttA 

O Hsieh Data 


• 5.0 -25 


0.0 2 S 5.0 

x/H 


75 10.0 


Test Cases 

Saiben Transonic Diffuser Flows 


/THROAT HEIGHT (H o ) = 44 mm 


INFLOW t.4 H 0 



81 

^ T'~ f 



g a sggi 




!■!•!■•] 










Static Pressure Distribution - N n Shock Case. R=0.862 


Top Wall 




Concluding Remarks 

• Developed MMLPG for adverse pressure gradient flows 

• Shown that MMLPG successfully computes boundary layer 
flows and transonic diffuser flows 

• Future work: continue validation 

modifications for separated flow 
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Turbulence Model Experiences for a 
Round-to-Rectangular Transition Duct 

Workshop on Computational 
Turbulence Modeling 

September 15-16, 1993 


Charles E. Towne 

NASA Lewis Research Center 
Cleveland, OH 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Outline 

• Proteus code description 

• Geometry and flow conditions 

• Numerical details 

• Turbulence model modifications 

• Convergence history 

• Computed results and 
comparison with experimental data 

• Concluding remarks 


TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Proteus Navier-Stokes Code 

• Reynolds-averaged, unsteady, compressible 
Navier-Stokes equations 

• Strong conservation-law form 

• Fully-coupled ADI solution procedure, 
Beam-Warming generalized time differencing 

• Second-order central differencing for 
spatial derivatives 

• Implicit steady/unsteady boundary conditions 

• Convection and diffusion terms linearized 
using second-order Taylor series expansion 

• Generalized nonorthogonal body-fitted coordinates 

• Baldwin-Lomax and Chien k-e turbulence models 


TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Proteus Navier-Stokes Code 

• 2-D, axisymmetric w/wo swirl, or 3-D flow 

• Thin-layer and Euler options 

• Wide variety of boundary conditions 

• Constant stagnation enthalpy option 

• Constant-coefficient or adaptive artificial viscosity 

• First- or second-order time differencing 

• Variety of time step selection methods 

• Output files for CONTOUR and PLOT3D 

• Highly vectorized for Cray computers 

• Extensively commented 

• Three-volume documentation set 


TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 



Circular-to-Rectangular Transition Duct 



STATION 1 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Flow Conditions 

• M ref = M inUt = 0. 2 to simulate 
incompressible experiment 

• Re ref = ( PURISM# = 195,000 

• Constant c p and c v , with y = 1.4 

• Constant at 519 °R 

• Constant stagnation enthalpy 


TRANSITION DUCT FLOW 
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Total Pressure Contours 




TRANSITION DUCT FLOW 


National Aeronautics and 
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Lewis Research Center 


Boundary and Initial Conditions 

• Boundary 

conditions 

Boundary 

Conditions Specified 

Inlet 

‘S 

II 

✓ N 

PO 

* 

K 

II 

© 

< 

li 

II 

O 

Exit 

m = map, u x = v x = =0 

Outer wall 

p r =0, 1/ = V = W = 0 

Centerline 

(p, u, v, w)cl = ave (p, u, v, w) CL+1 

0 = 0 

p r = u r = w r = 0, v = 0 

0 = 90 

p r = u r = v r = 0, w = 0 

• Initial conditions 

P = Poo. W = 

: (Uap)inlet. V = W = 0 


TRANSITION DUCT FLOW 
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Lewis Research Center 
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Numerics 

• 610) x31(0) x50(r) and 61 x31 x99 meshes 


• Local time step, with: 


Mesh 

Time Level 

CFL 

61x31x50 

1-1001 

1 


1001 - 2001 

10 


2001 - 3001 

10 


3001-4001 

20 

61x31 x 99 

1 - 1501 

1 


1501 -3001 

10 


3001 -4501 

10 


4501-6001 

15 


6001 -7501 

15 


• Cray Y-MP CPU time — just under 5 x 10 -5 
sec/grid point/time step, or 5.2 and 18.8 hours 
for the two cases 



TRANSITION DUCT FLOW 
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Computational Mesh 



Distance along coordinale line, s/R 


TRANSITION DUCT FLOW 
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Baldwin-Lomax Turbulence Model 

Inner region 

• ju t = pl 2 \Q\Re re f 
Outer region 

' Ml = KC cp pF KiebF wake Re ref 

where F wa ke = (y n)max^ max’ F max = rnax [j„l£2l(l — e )], 

and ( y n )max is the corresponding y n 

• Model based on boundary layer flows, without accounting 
for secondary flows and streamwise vorticity 

• At each streamwise station, the search region for was 
limited to the value of ( y n )max ’ n the y = 0 plane 

• Rationale — to prevent the secondary vortices, at the 
“ends” of the cross-section, from affecting the //, values 



TRANSITION DUCT FLOW 
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Effect of Turbulence Model Modification 







Distance along coordinate line, 


TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Effect of Turbulence Model Modification 




TRANSITION DUCT FLOW 
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Effect of Turbulence Model Modification 


Modified m mode! 


Original jj u mode! 








Computed mass (tux 


TRANSITION DUCT FLOW 
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Effect of Turbulence Model Modification 




TRANSITION DUCT FLOW 
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Mass Flow Rate Convergence 







TRANSITION DUCT FLOW 
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Convergence History 
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Peripheral Static Pressure Distribution 



- 1.0 - 0.6 - 0.2 0.2 0.6 1.0 
Circumferential distance, s/s re( 







Centerline Static Pressure Distribution 



Axia! distance, x/R 


TRANSITION DUCT FLOW 
Secondary Velocity Vectors 


National Aeronautics and 
Space Administration 
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TRANSITION DUCT FLOW 
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Space Administration 

Lewis Research Center 




TRANSITION DUCT FLOW 
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Total Pressure Contours 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 



Total Pressure Contours 




TRANSITION DUCT FLOW 
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Total Pressure Profiles Along Semi-Minor Axis 



Total pressure coefficient, (p r p;)/Q re t 




TRANSITION DUCT FLOW 
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Total Pressure Profiles Along Semi-Major Axis 



Total pressure coefficient, (p r pj)/Q ref 
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Concluding Remarks 

• The Proteus 3-D Navier-Stokes code has been used to compute 
the flow through a round-to-rectangular transition duct 

• The search region for F ^ in the Baldwin-Lomax turbulence model 
was limited to prevent excessively large //, values 

• The code correctly captures the basic physics — the generation 
of secondary flows and the resulting distortion 

• Agreement between computed and experimental results is good 
through the transition section 

• The computed secondary flows are overly damped in the 
downstream straight section 

• Further work is needed to investigate the effects of mesh 
resolution and turbulence model 
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Proteus experience with three different types of turbulence 

models 


Trong T. Bui 

NASA Lewis Research Center 
Cleveland, Ohio 


NASA Lewis Research Center INTERNAL FLUID MECHANICS DIVISION 

Cleveland, Ohio 


Contents 

1 . Overview of turbulence models used 

2. Practical issues 

3. Suggestions 


NASA Lewis Research Center 
Cleveland. Ohio 


NTERNAL FLUID MECHANICS DIVISION 


Overview of turbulence models used 


1 . Algebraic model 

• Baldwin-Lomax (BL) 

2. One-equation model 

• Baldwin-Barth one-equation (BB) 

3. Two-equation models 

• Chien k-s (CH) 

• Launder-Sharma k-e (LS) 


NASA Lewis Research Center 
Cleveland, Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Contents 

1 . Overview of turbulence models used 

2. Practical issues 

• CPU time requirement 

• Initialization of multi-equation models 

• y + dependency 

• Compressibility corrections 

• Where the turbulence models work well, and where they do not 
work too well 

3. Comments 


NASA Lewis Research Center 
Cleveland. Ohio 


INTERNAL fluid mechanics division 



uhu time requirement 


Not a major issue in using multi-equation turbulence models: 


{CPU time in sec/iter/grid point) 


Turbulence 

models 

Incomp, flat 
plate 

Comp, flat 
plate 

S-duct 

3-D shock/ 
B.L. inter. 

B.- L. 

6.031E-5 

9.042E-5 

5.168E-5 

7 . 888E-5 

B.- B. 

6. 650E-5 

9.436E-5 

N/A 

N/A 

Chien 

6.822E-5 

9.836E-5 

6 . 50 6E-5 

9.757E-5 

L.- S. 

8.667E-5 

11.77E-5 

N/A 

N/A 


• For flat plates, the Chien k-e model takes only about 10% more CPU 
time than the B.- L. model 

• For 3-D S-duct and shock/B.L. inter., the Chien k-e model takes about 
25% more CPU time than the B.- L. model 


NASA Lewis Research Center 
Cleveland, Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Initialization of multi-equation models 

Initial profiles of turbulence variables are needed to start the time 
marching process. The following steps have been found to work well: 

1 . |i t is computed using the BL model. 

2. BB model - k 2 /e is computed using m and the BB formula. 

3. k-e models - e is calculated from the local equilibrium assumption, 
and k is calculated using £, ji t , and the CH damping function f R . 


NASA Lewis Research Center 
Cleveland, Ohio 
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Initial k Profile for Glancing 
Shock Wave/TBL Interaction 



NASA Lewis Research Center INTERNAL FLUID MECHANICS DIVISION 

Cleveland, Ohio 


Y + dependency 


y + is needed to feed the Van Driest style damping functions in the near 
wall region. Not desirable for two main reasons: 

1 . Geometrical considerations: 

• Need to identify the nearest solid walls. This is difficult for complex 
3-D geometries with multiple solid boundaries. 

• Need to compute distances to solid walls. This is difficult to do for 
unstructured grids. 

• What is y + in base/wake flows? 

2. Ambiguity in the definition of y + for variable-property flows: 



Kv R 'r 

w 

\ Pi 


• Local or wall values for v-, and p 1 ? 

However, models that use y + appear to be more stable numerically than 
the ones that do not. 

INTERNAL FLUID MECHANICS DIVISION 


NASA Lewis Research Center 
Cleveland, Ohio 
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Effects of y + Computation, CH 




NASA Lewis Research Center INTERNAL FLUID MECHANICS DIVISION 

Cleveland. Ohio 
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Effects of Compressibility Corrections, CH 



0 5 000 10 e 1.000 10 7 1.500 10 7 2,000 10 7 2.500 10 7 3.000 10 7 

Re x 
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Contents 


1 . Overview of turbulence models used 


2. Practical issues 

• CPU time requirement 

• Initialization of multi-equation models 

• y + dependency 

• Compressibility corrections 

-► • Where the turbulence models work well, and where they do not 
work too well 

3. Comments 


NASA Lewis Research Center 
Cleveland, Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Where the turbulence models work well 


Attached, thin layer flows with little or no pressure gradients 






U/U, T/T # 


NASA Lewis Research Center 
Cleveland, Ohio 
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Where they don’t work too well 

Separated/Reattached flows and flows with large pressure gradients 





NASA Lewis Research Center INTERNAL FLUID MECHANICS DIVISION 

Cleveland, Ohio 


Comments 

• Two-equation k-s models can be used without serious CPU time 
penalty. However, no real gain in predictive capability compared with 
algebraic formulations. 

• To improve the predictive capability, need to extend the two-equation 
formulations with algebraic Reynolds stress models or anisotropic k-s 
models. 

• y + dependency limits the generality of turbulence models, and it should 
be removed if possible (preferably without significant cost in the 
stability and accuracy of the models). 

• Models are needed for the Reynolds heat flux terms. 


NASA Lewis Research Center 
Cleveland, Ohio 


INTERNAL FLUID MECHANICS DIVISION 
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Several Examples Where Turbulence Models 
Fail in Inlet Flow Field Analysis 


Bernhard H. Anderson 
NASA Lewis Research Center 
Cleveland, OH 


Inlet Flow Field Analysis 

Computational Uncertainties 

• Turbulence Modeling for 3D Inlet Flow Fields 

(1) Flows Approaching Separation 

(2) Strength of Secondary Flow Field 

(3) 3D Flow Predictions of Vortex Liftoff 

(4) Influence of Vortex-Boundary Layer Interactions 

• Vortex Generator Modeling 

(1) Representation of Generator Vorticity Field 

(2) Relationship Between Generator and Vorticity Field 


Inlet Flow Field Studies 

Goals and Objectives 


To advance the understanding, the prediction, and the control 
of intake distortion, and to study the basic interactions that 
influence this design problem. 

• To develop an understanding of and predictive capability 
for the aerodynamic properties of intake distortion and 
its management. 

• To establish a set of design guidelines to maximize the 
effectiveness of vortex flow control for the management, 
intake distortion 



Inlet Flow Field Benchmark Data Sets 

Turbulence and Vortex Generator Modeling 


• Fraser Flow A, Stanford Conference 1968 

• 727/JT8D-100 S-Duct Confirmation Experiment, 1973 

• Univ. Tennessee Diffusing S-Duct Experiments, 1986 & 1992 

• Univ. Washington TD410 Transition Duct Experiment, 1990 

• M2 129 Imperial College Laser-Doppler Experiment, 1990 

• DRA-Bedford Experiments on the M2 129 Intake S-Duct 

(2) DRA Surface Pressure and Engine Face. Experiment, 1989 

(3) DRA Phase IB Hot-Wire Flow Experiment, 1990 

(4) DRA Phase 2 Yawmeter Flow Experiment, 1991 

(5) DRA Phase 3 Vortex Generator Experiment, 1992 

• TD118 Bi-Furcated Transition Duct Experiment, 1994 


Reduced Navier-Stokes Analysis 

RNS3D Computer Code 

Velocity decomposition approach, Briley and McDonald (1979 & 1984) 

(1) Conservation form of the vorticity transport equation 

(2) Mass flow conservation, m = \j)u p dA = constant 
Non-orthogonal coordinate system, Levy, Briley and McDonald (1983) 
Arbitrary geometry gridfile description, Anderson (1990) 

(1) Recluster existing gridfile mesh distribution 

(2) Redefine the centerline space curve 

(3) Alter cross-sectional duct shape 
McDonald Camarata turbulence model 



Full Navier-Stokes Analysis 

PARC3D Computer Code 


• Originally developed by Pulliam & Steger as AIR3D (1980) 

(1) Conservation form of the governing equations 

(2) Beam & Warming approximate factorization algorithm 

(3) Central differencing within a curvilinear system 

• Addition of Jameson artificial dissipation by Pulliam, ARC3D (1981) 

• Developed for internal flow by Cooper, PARC3D (1987) 

(1) Baldwin-Lomax turbulence model 

(2) Diewert approximation to turbulence model in the 
reverse flow region of flow field 


Fraser Flow A, Stanford Conference 1968 

Geometry and Mesh Definition 



Distance from the Wall, inches 


Fraser Flow A, Stanford Conference 1968 

Comparison of Turbulence Models 
nz = 49, y + = 1.17 


RNS Analysis, MC Model 
FNS Analysis, BL Mode! 



Fraser Flow A, Stanford Conference 1968 
Comparison of turbulence Models 
nz = 49 ,y + = 1.17 


RNS Analysis, MC Model 
FNS Analysis, BL Model 
FNS Analysis, (k - e) Mod,eI 


O Data 


OO 


O -Q — &£ L. 


(a) Incompressible Shape Factor (HI) Distribution 



Axial Centerline Distance, ft. 

(b) Wall Skin Friction (Q) Distribution 



Fraser Flow A, Stanford Conference 1968 

Effect of y + on Flow Field Solution 



5 

4 

3 

2 

1 

0 


— nz = 49, y+ = 0.584 

— nz ~ 49, = 1.168 

nz = 49, y+ = 5.739 

O Data 


-Q .0—0*- 


oo 


(a) Incompressible Shape Factor (HI) Distribution 



(b) Wall Skin Friction (Q) Distribution 


RNS Analysis, McDonald-Camarata Model 


Fraser Flow A, Stanford Conference 1968 

Effect of Mesh Resolution on Flow Field Solution 


- nz = 199, y + = 0.477 

- nz = 99, y* = 0.520 

- nz = 49, y * = 0.584 
O Data 



(a) Incompressible Shape Factor (HI) Distribution 



(b) Wall Skin Friction (Q) Distribution 

RNS Analysis, McDonald-Camarata Model 



Fraser Flow A, Stanford Conference 1968 

Conclusions 

(1) The current generation of turbulence models were unable to 
predict the complete state of the diffuser boundary layer 
approaching flow separation. 

(2) Both near wall and mesh resolution separately played an 
important role in accurate solutions to wall skin friction 
distribution in flows characterized as "approaching separation", 
but had little effect on the solution for the incompressible 
shape factor development. 

(3) It is important that grid independent solution be demonstrated 
before judgements about the turbulence models be stated, and 
that the complete state of the wall boundary layer be 
considered within this evaluation. 


727/JT8D-100 Inlet S-Duct 

Geometry and Mesh Definition 



727/JT8D-100 Inlet S-Duct 

Engine Face Flow Field 
Generator Config. 12 



Experiment Analysis 

With Engine Dome Without Engine Dome 


RNS Analysis, McDonald-Camarata Model 


727/JT8D-100 Inlet S-Duct 

Engine Face Ring Distortion Characteristics 
Generator Config. 12 


O Experiment (With Engine Dome) 
_ Analysis (Without Engine Dome) 



Engine Face Radius, Inches (Model Scale) Engine Face Radius, Inches (Model Scale) 


Radial Distortion 60° Sector Circumferential Distortion 


RNS Analysis. McDonald-Camarata Model 



727/JT8D-100 Inlet S-Duct 

Conclusions 


(1) The current turbulence models predict the overall performance 
level of vortex generator installation remarkably well, although 
much of the detailed flow structure was not resolved. 

(2) Turbulence models in 3D inlet flow field analysis can also be 
evaluated on the basis of standard engine performance parameters, 
such as radial and circumferential ring distortion descriptors, 
which provide a sensitive discriminator measuring the state of 

the overall compressor face flow field. 


Univ. Tennessee Diffusing S-Duct 

Geometry and Mesh Definition 





Univ. Tennessee Diffusing S-Duct 

Topology of Vortex Liftoff 



Experiment 


RNS Analysis. McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 

Topography of Vortex Liftoff 



RNS Analysis, McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 
Total Pressure Coefficient Contours 
Without Vortex Generators 




Analysis 


RNS Analysis, McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 
Total Pressure Coefficient Contours 
With Vortex Generators 



Experiment Analysis 


RNS Analysis, McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 

Effect of y + on Engine Face Flow Field 
Total Pressure Coefficient Contours 



RNS Analysis. McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 

Effect of y + on Circumferential Distortion 


y + = 0.5 
y + — 8.5 



RNS Analysis, McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 

Conclusions 

(1) Initial value space marching 3D RNS procedures adequately 
described the topological and topographical features of 3D flow 
separation associated with vortex liftoff. 

(2) The current turbulence models predicted the overall structure of 
vortex generator installation remarkably well, although much of the 
detailed flow structure was not resolved. 

(3) Adequate near wall resolution was necessary to obtain an accurate 
solution of the phenomena of vortex lift-off. 

(4) Circumferential ring distortion is a sensitive discriminator in 
measuring the state of the engine face flow field. 


Univ. Washington TD410 Transition Duct 

Geometry and Mesh Definition 


Univ. Washington TD410 Transition Duct 

Case Definitions 


Case 

Grid 

Total 

CPU (min) 

7+ 

td410.$ 

199 x 121 x 521 

12,545,159 

149.0 

0.565 

td410.1 

99 x 121 x 521 

6,241,059 

73.4 

0.623 

td410.3 

49 x 121 x 521 

2,851,849 

36.0 

0.568 

td411.1 

99 x 91 x 521 

4,693,689 

57.1 

0.623 

td412.$ 

199 x 61 x 521 

6,324,419 

72.9 

0.565 

td412.1 

99 x 61 x 521 

3,146,319 

38.2 

0.623 

td412.3 

49 x 61 x 521 

1,557,269 

18.7 

0.568 

td412.4 

49 x 61 x 401 

1,198,589 

14.7 

0.568 


Univ. Washington TD410 Transition Duct 

Effect of Radial Grid Resolution 



(a) Velocity Profiles Along the Minor Axis Surface 



(b) Velocity Profiles Along the Major Axis Surface 


RNS Analysis, McDonald-Camarata Model 







Univ. Washington TD410 Transition Duct 

Effect of Radial Grid Resolution 
Wall Skin Friction 



(a) Skin Friction C/ along the Minor Axis Wall 



Axial Centerline Distance, X/Ri 

(b) Skin Friction Q along the Major Axis Wall 


RNS Analysis, McDonald-Camarata Model 


Univ. Washington TD410 Transition Duct 

Comparison with Experimental Data 
Velocity Contours 



RNS Analysis, McDonald-Camarata Model 



Univ. Washington TD410 Transition Duct 

Conclusions 

(1) The current generation of turbulence models predicted the 
overall development of vortex formation reasonably well, 
although there were important discrepencies which could not 
not be explained as inadequate near wall or mesh resolution. 

(2) Radial mesh resolution had the largest impact in the region 

along the major axis where the vortex pair was formed. 

(3) It is important that grid independent solution be demonstrated 
before judgements about the turbulence models be stated. 

(4) Fully 3D grid independent solutions were achieved with a 
Reduced Navier-Stokes analysis. 


DRA M2129 Diffusing Inlet S-Duct 

Geometry and Mesh Definition 


DRA M2129 Diffusing Inlet S-Duct 

Performance Characteristics 
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- RNS Analysis, MC Model 
• FNS Analysis, BL Model 
O Data 


Engine Face Total Pressure Recovery 



Q.2 0A~ (L6 (T8 Flo 


Throat Mach Number, M t 
(b) Engine Face DQo Distortion 


DRA M2 129 Diffusing Inlet S-Duct 

Separation Characteristics 


— RNS Analysis, MC Model 
• FNS Analysis, BL Model 



Throat Mach Number, M t 



Static Pressure Ratio, pjPu 


DRA M2129 Inlet S-Duct 

Wall Static Pressure Distribution 
AGARD Test Case 3.1 


- RNS Analysis, MC Model 

- FNS Analysis, BL Model 



.6 - 
.5 - 

A L 

(a) 0 = 0.0° Surface Element. 



(b) 0 = 180.0° Surface Element. 


DRA M2129 Diffusing Inlet S-Duct 

Engine Face Distortion Characteristics 
AGARD Test Case 3.1 




(b) 60°-Sector Circumferential Distortion 


DRA M2129 Inlet S-Duct 

Conclusions 


(1) Both Full Navier-Stokes (FNS) and Reduce Navier-Stokes (RNS) 
analyses adequately describe the overall flow physics of vortex 
liftoff, but consistently predict the location of liftoff further 
downstream in the duct inlet than was indicated by data. 

(2) The current generation of turbulence models were unable 
to desribed the influence of separation on the main pressure 
field for "strong" vortex liftoff interactions. 

(3) The current generation of mixing length turbulence models 
give remarkable good performance results, while the existing 
discrepancies between data and analysis can be attributed 
primarily to the over prediction of the liftoff location. 


TD118 Bi-Furcated Transition Duct 

Geometry and Mesh Definition 


TD118 Bi-Furcated Transition S-Duct 

Research Objectives 


• To demonstrate diffuser duct technology advancement by using 
CFD to design a "conventionally shorter 7 ' transitioning S-duct 
configuration for application towards high speed inlet systems. 

• To develop a computational protocol whereby turbulence model 
evaluations can be made between different computer codes. 

• To develop a benchmark data set to evaluate CFD analysis and 
turbulence models, which cover fundamental flow phenomena as 
well as overall flow field physics as determined by standard 
engine performance parameters. 


TD118 Bi-Furcated Transition S-Duct 

Effect of Turbulence Model on Inlet Performance 


Analysis 

Turbulence Model 

PtavjPt o 

DH 

DCs o 

RNS3D 

McDonald-Camarata 

0.955 

0.167 

0.116 

PARC3D 

Baldwin-Lomax 

0.959 

0.163 

0.135 

PARC3D 

P.D. Thomas 

0.970 

0.153 

0.073 

PARC3D 

Launder-Spaulding 

0.975 

0.150 

0.177 


DH — (P/max Ptm\n)l Ptavc 


DCs 0 — {Ptave 




TD118 Bi-Furcated Transition S-Duct 

Effect of Turbulence Model on Engine Face Flow Field 



(a) RN$ Analysis, McDonald-Camarata Model 


(c) FNS Analysis, P. D. Thomas Model 



(b) FNS Analysis, Baldwin-Lomax Model (d) FNS Analysis, Launder-Spalding (/< — £) Model 


TD118 Bi-Furcated Transition S-Duct 

Effect of Turbulence Model on Wall Boundary Layer 
6 = 90.0° Surface Element 


Analysis 

FNS, BL Model 

FNS, PDT Model 




Axial Centerline Distance, XJD < 


(b) Wall Skin Friction C { Distribution. 


TD118 Bi-Furcated Transition Duct 

Effect of Turbulence Model on Ring Distortion 


Analysis 

— FNS, BL Model 




Inlet Flow Field Analysis 

Concluding Remarks 

(1) Difficulties in complex 3D flow fields often arise because 
fundamental 2D aerodynamic interactions have not been 
adequately resolved. 

(2) Near wall (y + ) and radial mesh resolution ( nz ) play an 
important role in fundamental 2D and complex 3D flow 
field analysis. 

(3) Judgements about turbulent models should not be stated until 
grid independent solutions have been established. 

(4) Adequateness of turbulence models in inlet flow field analysis 
should also be made on the basis of fundamental performance 
parameters used to quantify the "goodness" of the flow 
entering the engine. 



